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SECTION  1 


INTRODUCTION 

Kaman  Sciences  Corporation  (KSC) ,  under  contract  to  the 
Defense  Nuclear  Agency  (DNA) ,  has  reviewed  the  methodology 
developed  by  United  Kingdom  (UK)  Ministry  of  Defence  and  used  by 
the  U.S.  Army  for  determining  the  thermal  transmissivities  in 
Northwest  Europe  for  tactical  nuclear  weapons.  Kaman  then 
compared  the  UK  methodology  with  the  present  KSC  methodology 
which  is  being  employed  in  the  new  EM-I  (Capabilities  of  Nuclear 
Weapons)  for  Northwest  Europe.  Based  on  this  comparison,  Kaman 
Sciences  Corporation  recommended  changing  to  a  methodology  which 
would  update  the  present  data  base  for  determining  the  thermal 
transmissivities  for  tactical  nuclear  weapons  employed  in 
Northwest  Europe. 

1.1  OUTLINE  OF  REMAINING  CHAPTERS. 

Chapter  2  presents  a  discussion  of  the  procedures  employed 
by  Maj.  B.M.  Cooke  in  generating  thermal  transmissivity  data  for 
the  UK  Ministry  of  Defence.  The  chapter  includes  a  discussion  of 
the  basic  transmission  equation,  the  meteorological  data, 
geographical  data,  and  data  interpolation  and  extrapolation  used 
in  the  UK  methodology.  Chapter  3  offers  a  presentation  of  the 
KSC  methodology  employed  for  this  effort.  This  section  includes 
a  description  of  the  RECIPE  code,  the  transmission  data  base, 
weather  analysis  requirements,  transmission  techniques,  and 
probabilistic  transmission  employed  by  KSC.  Chapter  4  presents  a 
comparison  of  the  results  using  the  UK  versus  the  KSC  methodo- 
lociy.  Chapter  5  describes  the  QRT  data  base.  This  section  also 
discusses  the  limitations  in  the  QRT  calculational  approach  and 
data  base.  Chapter  6  offers  a  recommended  approach  for  revisions 
in  the  QRT  data  base  for  use  in  the  KSC  methodology. 
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SECTION  2 


THE  UK  APPROACH 

The  energy  flux  of  thermal  radiation  from  a  nuclear  fireball 
decreases  rapidly  with  increasing  distance  from  the  burst.  This 
is  due  not  **ly  to  inverse-distance-squared  geometric  attenua¬ 
tion,  but  also  to  the  absorption  and  scattering  of  the  radiation 
during  its  passage  through  large  distances  in  air.  This  latter 
effect  is  termed  atmospheric  transmission.  In  the  wavelength 
region  of  interest,  0.3  to  4.0  ji  ,  atmospheric  transmission 
depends  on  such  factors  as  the  absorption  of  infrared  radiation 
by  water  vapor  and  carbon  dioxide,  the  absorption  of  ultraviolet 
and  visible  radiation  by  ozone,  multiple  scattering  at  all 
wavelengths  from  solid  particles  suspended  in  air  and  air  mole¬ 
cules,  reflections  from  the  ground  and  clouds,  etc.  All  of  these 
factors  vary  with  altitude  and  with  prevailing  weather  condi¬ 
tions.  Interactions  such  as  scattering,  which  is  produced  by  the 
reflection  and  refraction  of  light  rays  by  certain  atmospheric 
constituents  such  as  dust,  smoke,  and  fog,  divert  the  rays  from 
their  original  paths  and  result  in  diffuse,  rather  than  direct 
transmission  of  the  radiation.  Therefore,  a  receiver  with  a 
large  field  of  view  (that  is,  most  military  targets)  receives 
radiation  that  has  been  scattered  toward  it  from  many  directions 
in  addition  to  the  directly  transmitted  radiation.  Because  the 
mechanisms  of  absorption  and  scattering  depend  on  the  wavelength, 
the  atmospheric  transmittance  depends  not  only  upon  the  atmos¬ 
pheric  conditions,  but  also  upon  the  spectral  distribution  of  the 
radiation  emitted  by  the  fireball. 

2.1  THE  RILEY  EQUATION. 

The  approach  employed  by  the  UK  was  based  on  an  empirical 
equation  for  thermal  transmission  in  the  atmosphere  developed  by 
George  Riley  (Reference  1).  In  1970,  Maj.  C.B.  Pritchett 
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continued  work  originally  started  by  B.  Rankin  on  the  prediction 
of  thermal  transmissivity  factors  for  Northwest  Europe,  using 
meteorological  data  compiled  over  a  period  of  some  10  years  at 
airfields  in  Northwest  Germany.  Riley's  formula  was  used  to 
calculate  transmissivity  values  which  would  not  be  exceeded  at 
least  95,  50,  and  5  percent  of  the  time.  Maj.  Cooke  (Reference 
2)  extended  this  work  and  offers  the  following  rationale  for  the 
selection  of  this  methodology. 

"Several  relationships  are  available  for  the  prediction  of 
thermal  transmission  factors  from  known  visibility  conditions. 
Eldridge  (Reference  3)  presented  some  of  these,  one  of  which  was 
Riley's  formula.  Riley  has  derived  functional  relationships 
between  thermal  transmissivity  and  the  distance  between  the 
source  and  the  receiver  for  various  meteorological  conditions. 
This  method  was  selected  because: 

a.  It  is  applicable  to  the  climate  of  Northwest  Germany  as 
it  includes  the  effects  of  cloud  cover. 

b.  The  ranges  (0.2  to  25  km)  over  which  the  relationship 
can  be  used  are  suitable  for  military  applications. 

c.  The  visibility  range  acceptable  to  the  relationships 
(1.3  to  186  Km)  is  sufficiently  wide  to  encompass  the  highest  and 
lowest  probabilities  likely  to  be  used. 

d.  They  are  derived  for  surfaces  normal  to  the  thermal 
fluence,  which  is  assumed  to  be  relevant  to  military  operations." 

Riley's  relationship  allows  transmissivity  to  be  predicted 
from  cloud  base  height,  visibility,  and  snow  cover  data.  Unfor¬ 
tunately,  Cooke  was  unable  to  obtain  data  which  interrelates 
these  three  variables.  It  was  only  possible  to  obtain  data 
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relating  the  frequency  of  cloud  base  height  to  visibility. 
Therefore,  the  effects  of  snow  cover  were  ignored. 

Riley  at  the  Air  Force  Cambridge  Research  Laboratory  has 
developed  an  empirical  transmission  prediction  technique  by 
representing  the  transmission  of  visible  radiation  as: 

T  -  e  ”0D  +  D2(Ae  _3*lD  +  Be  -1*12D  +  ce  “  aD)  (1) 

whero  T  if  the  transmission  over  the  path  length  D  (km)  ,  0  is  an 
attenuation  coefficient,  and  A,  B,  C  and  CC  are  parameters 
depending  on  the  meteorological  conditions.  The  parameter  values 
presented  are  based  on  data  taken  over  path  lengths  from  0.2  to 
25  km  with  values  of  0 D  between  0.02  and  40.  The  experimental 
data  were  obtained  from  References  5,  6,  7,  8,  and  9.  Riley 
suggested  determining  the  value  of  0  wherever  possible  by 
measuring  the  direct  radiation.  The  alternative  used  by  Cooke 
is: 


0  -  3.91/V  (2) 

where  v  is  the  visual  range  defined  as  the  range  where  the 


transmission  decreases  to  0.02. 
is  the  meteorological  range. 

The  values  of  the  variables 
equation  2.1  are: 


The  correct  terminology  for  p 


as  applied  by  Maj.  Cooke  in 


1.29  -  for  an  overcast  sky  less  than  1500  meters  cloud 

B  *  base  height, 

0.151  -  clear  sky. 
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c 


f 0.0102  -  for  an  overcast  sky  below  300  meters  cloud 
base  height, 

-  0.033  -  for  an  overcast  sky  with  a  cloud  base  height 

between  300  and  1500  meters, 

'0.0033  -  clear  sky. 


a  « 


fO.14500*045  -  clear  sky, 

0  12 

0 . 220w* -  for  an  overcast  sky, 


where  an  overcast  sky  has  greater  than  five-eighths  cloud  cover. 


This  expression  applies  to  horizontal  paths  near  sea  level 
and  for  a  20  detector  aimed  at  the  source. 


The  first  terms  in  this  expression  represent  the  direct 
radiation  while  the  other  terms  account  for  scattered  and 
reflected  radiation. 

2.2  THE  METEOROLOGICAL  DATA. 


Cooke  was  supplied  visibility  and  cloud  base  height  data 
taken  at  intervals  of  three  hours  over  a  period  of  approximately 
ten  years  at  Gutersloh,  Fassberg,  Hannover,  Hamburg,  and 
Dusseldorf.  Annual  and  monthly  means  were  also  provided  and 
formed  the  basis  of  the  UK  analysis.  Annual  and  monthly  statis¬ 
tics  were  assembled  which  related  the  frequency  of  the  three 
categories  of  cloud  base  height  and  thirteen  categories  of 
visibility.  The  three  categories  of  cloud  base  heights  are: 

Group  1,  less  than  300  meters;  Group  2,  300  to  1,500  meters;  and 
Group  3,  above  1,500  meters.  Cooke  combined  the  case  of  no 
clouds  into  Group  3.  The  thirteen  categories  of  visibility  were 
selected  ranges  from  0  to  200  meters  up  to  8,000  to  16,000  meters 
as  shown  in  Table  1.  Tables  were  then  drawn  up  of  the  number 
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Table  1.  Data  Collection  Format  for  Visibility /Cloud  Base  Height  Data. 
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of  occurrences  of  various  weather  conditions  for  the  set  of 
annual  means  and  for  each  set  of  monthly  means.  From  these,  the 
cloud-base  height /visibility  distributions  were  derived. 

2.3  GEOGRAPHICAL  VARIATIONS. 

It  was  not  known  how  dependent  the  weather  statistics  were 
on  the  geographical  variations  of  the  recording  stations.  It 
would  be  simpler  if  a  single  set  of  statistics  could  be  used  to 
represent  the  cloud  base/visibility  relaticnships  of  all  five 
locations.  However,  if  there  were  significant  differences 
between  the  weather  statistic?  of  the  five  locations,  then 
reservations  would  need  to  be  placed  on  the  use  of  a  single  set 
of  weather  statistics  to  represent  Northwest  Germany. 

A  chi-squared  value  was  used  by  Cooke  as  a  measure  of 
association  for  the  data  containing  the  variables  location  and 
visibility.  The  larger  the  chi-squared  value,  the  greater  is  the 
independence,  and  conversely,  the  smaller  the  chi-squared  value, 
the  greater  the  degree  of  association. 

In  order  to  determine  the  number  of  occurrences  of  various 
weather  conditions  at  the  various  geographical  locations  from  the 
meteorological  data  over  the  period  of  ten  years,  tables  were 
drawn  up  for  each  set  of  annual  mean  and  monthly  mean  weather 
data  and  for  each  cloud  base  height  group,  as  shown  in  Table  2. 
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Table  2.  Data  Format  for  Chi-Squared  Test. 


LOCATION  VISIBILITY 

CATEGORY  CATEGORY 

"S"  "d" 


Dusseldorf 

1 

nllij  n12ij  n13ij 

nldij 

Fossberg 

2 

n21ij  n22i j 

n2dij 

Guterslok 

3 

n31ij 

n3dij 

Hamburg 

4 

n41i  j 

n4dij 

Hannover 

5 

n51ij 

n5dij 

Where  ngdij  is  the  total  number  of  occurrences  of  the  given 
weather  conditions  over  the  ten  year  period,  the  suffixes  have 
the  following  significance:  "S"  is  the  visibility  category  coded 
into  groups  1  to  13,  where  1  is  the  visibility  of  column  A  of 
the  meteorological  data  (Table  1),  i.e.,  less  than  200  meters;  2 
is  the  visibility  of  column  B  (200  to  400  meters),  etc.,  up  to 
13,  corresponding  to  column  M  of  the  weather  data  (8,000  to 
16,000  meters) . 

"i"  is  the  cloud  base  height  group  coded  into  the  three 
categories: 


1,  less  than  300  meters 

2,  between  300  and  1,500  meters,  and 

3,  above  1,500  meters. 

"j"  refers  to  the  particular  set  of  weather  data  of  which  there 
are  13  groups: 
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1 ,  the  annual  means 

2,  the  means  for  January, 

3,  the  means  for  February,  etc.,  up  to 
13,  the  means  for  December. 


Two  specimen  tables  are  reproduced.  Table  3  gives  the 

number  of  occurrences  n_ . . .  for  the  annual  means  and  Table  4 

Sai] 

gives  the  number  of  occurrences  for  the  January  means.  For  a 
table  of  "S"  rows  and  "d"  columns,  the  number  of  degrees  of 
freedom  for  chi-squared  is  (S-l)(d-l).  For  each  set  of  weather 
statistics,  ij  chi-squared  is  calculated.  The  results  for  the 
annual  means  and  the  January  means  are  also  shown  in  Tables  3  and 

4.  The  complete  set  of  chi-squared  values  for  all  the  months  is 
given  in  Table  5.  For  the  means  annual  weather  statistics,  the 
chi-squared  values  show  a  high  degree  of  association  between  the 
five  locations.  Because  of  this  homogeneity,  it  is  further 
assumed  that  the  means  statistics  can  represent  the  annual 
weather  data  for  these  five  locations.  Moreover,  as  the 
locations  are  widely  spread,  it  is  considered  that  one  set  of 
weather  statistics  can  reasonably  describe  the  annual  weather 
conditions  in  Northwest  Germany  in  areas  which  are  similar  to  the 
geophysical  conditions  of  those  five  locations,  i.e.,  fairly  low 
lying  and  flat. 


For  the  analysis  of  the  monthly  statistics,  for  cloud  base 
heights  Groups  2  (between  300  and  1,500  meters)  and  3  (above 
1,500  meters),  there  was  again  a  high  degree  of  association  for 
the  five  locations.  For  cloud  base  height  Group  1  (less  than  300 
meters),  there  was,  however,  a  diminution  of  the  evidence  for 
association,  particularly  for  the  months  April  to  October. 

Indeed,  the  value  of  chi-squared  for  August  indicates  some 
positive  degree  of  independence  between  the  data  provided  by  the 
various  weather  stations. 
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Table  3.  Annual  Statistics  (j-1)  Number  of  Occurrences  nsdjj. 
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Table  4.  January  Statistics  (j-2)  Number  of  Occurrences  n8dij. 
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Table  5.  Results  of  the  Tests  of  Geographical  Association. 


Parted 

<J> 

Cloud-late 

Height 

Croup 

(1) 

Chi-Squared 
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N usher 
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Oegreea 
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li.a 

1 

Annual 

2 

».o 

48 
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1 
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mm 
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48 
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10 
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1 

I 
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17.2 

48 

17 

3 

13.7 

18 

1 

77.7 
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2 

19.4 

48 
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3 
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tmm 

October 

2 

21.3 

48 

"El' 

3 

19.3 

KB 

1 

29.0 

34 

November 

2 

24.1 

48 

33 

3 

19.4 

34 

1 

u 

37 

December 

2 

48 

3 

MSM 

_!LJ 

NOTES: 

1.  The  values  of  chi— squared  for  the  mean  annual  statistics 
are  all  within  the  99.5%  level. 

2.  The  values  of  chi-squared  for  the  monthly  statistics  for 
cloud-base  height  Groups  2  &  3  are  all  within  the  99%  level. 

3.  The  values  of  chi-squared  for  the  months  April  to  October 
for  cloud-base  height  Group  1  were  not  significant  at  the  75% 
level. 

4.  The  value  of  chi-squared  for  the  month  of  August  for  cloud 
base  height  Group  1  was  not  significant  at  the  0.5%  level. 
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Some  nil  reports  occurred  in  a  number  of  the  monthly  weather 
statistics.  Further  tests  were  conducted  to  investigate  the 
possibility  of  these  causing  errors  in  the  statistical  analysis 
of  association.  The  weather  data  were  rearranged  into  a 
different  set  of  visibility  categories  to  eliminate  nil  reports 
and  the  chi-squared  values  were  recomputed.  The  results  obtained 
by  Cooke  did  not  differ  appreciably  from  the  original  results 
present  previously. 

2.4  DATA  INTERPOLATION  AND  EXTRAPOLATION. 


The  monthly  and  annual  weather  data  available  to  describe 
cloud  base  height  and  visibility  were  far  from  complete.  In 
Figure  2.1  the  known  mean  annual  distributions  of  the  three  cloud 
base  heights  chosen  and  the  overall  mean  are  shown.  It  can  be 
seen  that  little  long  range  visibility  data  were  available  for 
cloud  base  heights  Groups  2  and  3.  This  information  was  needed 
to  predict  the  various  probability  levels  of  transmissivity. 


The  annual  visibility  curves  in  Figure  1  appear  to  be  of 
cumulative  Gaussian  form.  For  each  curve  there  is  suggested  a 
relationship  of  the  standardized  form 


P 


.(x-  Mij)/  0ij 


1 

”1? 


exp(-*jt*)dt 


(3) 


where 


P  *  cumulative  occurrence, 

X  *  visibility  measured  on  a  logarithmic  scale, 
j  «  mean  vi  sibility  measured  on  a  logarithmic  scale, 
o_  ■  standard  deviation, 


and,  as  before. 
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Figure  1.  Measured  Cloud  Base  Height  Visibility  and  Mean 
Visibility  Distributions,  At.nual 


i  =  1  for  cloud  base  height  Group  1,  etc., 
j  *  1  for  the  annual  mean  weather  data, 
j  *  2  for  January  Mean  weather  data,  et:c. 

The  cumulative  Gaussian  curve  would  normally  be  expected  to 
operate  from  0  to  100  percent  of  recorded  data.  However,  it  can 
be  seen  from  Figure  1  that  cloud  base  height  Group  1  appears  to 
flatten  out  at  around  7  percent.  Furthermore,  an  analysis  of  all 
the  annual  and  monthly  statistics  showed  that  the  minimum 
probability  levels  for  cloud  base  Group  1  tailed  off  between  5 
and  2)  percent,  and  for  cloud  base  Group  2  and  3,  between  0  and  '' 
percent. 


It  was  assumed  by  Cooke  that  this  early  flattening  out  is 
produced  by  factors  which  solely  affect  the  low  visibility 
levels.  One  such  factor  for  the  cloud  base  height  Group  1  could 
arise  when  the  cloud  base  height  reaches  ground  level.  For  the 
higher  cloud  base  height  Groups  2  and  3,  industrial  haze  may 
introduce  factors  affecting  low  visibility  readings.  Moreover, 
the  general  practice  of  recording  the  worst  visibility  in  any 
direction  is  likely  to  increase  the  incidence  of  recorded  low 
visibility  readings  and  it  suggests  that  there  is  a  marked 
probability  of  these  being  nil  visibility.  In  such  circum¬ 
stances,  it  becomes  unrealistic  to  establish  a  Gaussian  fit  from 
0  to  100  percent  incidence  of  visibility  reading  and  the  above 
equation  is  rewritten  to  establish  a  model  of  the  form 


P-M . 
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1 

27 


J*.  •)/ 
13 


a.  . 

13 


(-»it2)dt 


(4) 


where  ^ 


optimum  probability  le/el  of  recorded  visibilities. 
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By  using  probit  analysis,  the  equation  can  be  further  re¬ 
written 


Gij  x  +  Aij- 


(5) 


where  Y  is  the  probit  corresponding  to  P-M^j 

with  »  1/  o  and,  in  order  to  avoid  the  necessary  computing 
with  negative  numbers,  the  customary  step  is  taken 


Au 


K  - 


(6) 


where  K  ■  5. 


a  .  ■ 


Computer  programs  where  used  to  obtain  the  values  of  G^j  and 
Aij  for  the  annual  and  monthly  cloud  base  height/visibility 
distributions. 


(1)  An  iterative  method  was  used  to  obtain  the  maximum 
likelihood  estimates  for  the  parameters  A^  j  and  G^  in  equation 
(2.5)  for  each  set  of  statistics. 


(2)  The  iteration  was  repeated  with  M^j  *  M^j  ±  5M 
until  there  is  the  maximum  agreement  between  observed  data  and 
data  predicted  by  equation  (4) .  A  chi-squared  test  was  used  as 
an  indicator  for  the  goodness  of  fit  of  the  cumulative  Gaussian 
curve  model. 


Samples  of  the  pictorial  comparisons  between  predicted  and 
observed  data  for  annual  and  monthly  cloud  base  height/visibili¬ 
ty  distributions  are  given  in  Figures  2  to  7  for  the  annual  data 
(for  all  three  cj. -  ad  base  height  groups)  .  A  reasonable  level  of 
agreement  between  the  cumulative  Gaussian  models  and  the  observed 
data  was  indicated  by  the  chi-squared  tests. 


15 


®ou#jjnooo  ®0»iu»3J8d  ®Aji®|numo 


Figure  2.  Visibility  Distribution  -  Cloud  Base  Height  Group  1. 
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Figure  3.  Visibility  Distribution  -  Cloud  Base  Height  Group  2 
NW  Germany  -  Annual. 
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Figure  4.  Visibility  Distribution  -  Cloud  Base  Height  Group  3 
NW  Germany  -  Annual. 
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Figure  5.  Visibility  Distribution  -  Cloud  Base  Height  Group  1  - 
NW  Germany  -  January. 


19 


100  1000  10000  100000 

*  Visibility.  Meties 


Figure 


100 


Figure  7.  Visibility  Distribution  -  Cloud  Base  Height  Group  3  - 
NW  Germany  -  January. 
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2.5  THERMAL  TRANSMISSION  FACTORS. 


From  each  Gaussian  model  of  the  annual  and  monthly  weather 
statistics,  an  entire  range  of  weather  data  was  predicted  and 
applied  to  Riley's  equation  to  obtain  thermal  transmissivities. 

Thermal  transmission  factors  were  computed  over  a  range  of 
0.25  to  23.75  km  at  0.5  km  intervals  for  various  probability 
levels,  and  they  are  presented  in  graphical  form  as  a  relation¬ 
ship  between  transmission  factors  and  distance  from  source  for  a 
series  of  probabilities.  The  mean  annual  transmission  factors 
are  given  in  Figure  8.  Similar  figures  were  computed  for  each  of 
the  monthly  mean  weather  statistics. 

The  relationship  with  the  lowest  transmissivities  represents 
the  5  percent  (i.e.,  on  5  percent  of  occasions  these  will  be  the 
transmissivity  described  or  less) .  The  next  highest  relationship 
is  10  percent  case  and  additional  relationships  are  given  in 
steps  of  10  percent  up  to  90  percent  case.  They  are  completed  at 
the  95  percent  case  (i.e.,  when  there  is  a  95  percent  proba¬ 
bility  that  the  given  transmission  factor  will  not  be  exceeded) . 

At  present,  thermal  calculations  for  military  equipment 
survivability  and  effects  tables  use  the  95,  50,  and  5  percent 
transmissivity  levels.  Comparing  the  mean  monthly  with  the  mean 
annual  transmissivities  for  the  three  probabilities  will  show 
some  seasonal  variations.  These  seasonal  variations  are  suffi¬ 
ciently  small  (10  percent)  to  be  ignored  in  the  95  percent 
probability  case,  but  the  maximum  variation  of  +75  percent  for 
the  5  percent  case  may  be  sufficiently  high  for  there  to  be  a 
need  to  recognize  seasonal  variations  in  certain  circumstances. 
However,  in  both  equipment  survivability  and  effects  tables 
studies,  for  the  poor  transmissivity,  it  is  the  worst  case  that 
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Figure  8.  Range  of  Probabilities  -  Annual 
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Distance,  km 


is  of  interest.  Here  the  variation  is  only  20  percent,  which  was 
viewed  as  insignificant  when  the  low  transmission  factor  (0.2) 
was  considered. 

2.6  CRITIQUE  OF  RILEY'S  EQUATION. 

In  the  analysis  performed  by  Kaman  Sciences  Corporation, 
several  concerns  were  raised  on  the  use  of  the  Riley  equation  for 
calculation  of  the  thermal  transmissivity. 

The  Riley  equation  is  an  empirical  fit  and  as  such  is  crude 
compared  to  presently  available  techniques.  The  relation  is 
based  in  part  on  a  series  of  experiments  using  a  tungsten  lamp. 
Therefore,  the  equation  can  not  be  affected  by  changes  in  the 
source  spectrum  which  would  result  from  the  use  of  different 
yield  weapons.  In  addition,  Riley  used  a  standard  surface  and 
does  not  account  for  various  terrain  features  (e.g.,  desert, 
forest,  or  urban  area)  which  may  apply  in  different  tactical 
scenarios.  The  relationship  is  valid  for  visible  light  and  does 
not  account  for  radiation  outside  the  visible  spectrum.  The  path 
radiation  travels  according  to  Riley's  equation  is  assumed  to  be 
horizontal  at  sea  level  with  a  constant  air  density.  This 
ignores  the  effect  of  varying  air  density  as  radiation  travels 
from  the  fireball  to  the  target  at  some  other  altitude,  which 
could  be  of  significance  for  both  land  and  air  targets.  The 
relationship  does  nov*  adequately  model  multiple  scattering 
contributions . 

There  also  exists  some  uncertainty  in  the  definition  in  the 
visibility  as  defined  by  Riley  and  employed  by  Cooke.  Riley 
defines  visual  range  applicable  to  his  relationship  as  the  range 
at  which  the  direct  beam  transmission  is  2  percent.  However, 
this  was  not  the  definition  employed  by  the  meteorological 
stations  which  used  an  "eyeball"  method.  Cooke  estimated  that 
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visibility,  commonly  observed  by  the  "eyeball"  method,  is  about 
one-half  of  the  visual  range.  Cooke  used  an  estimate  of  the 
relationship  between  visual  range  and  visibility  that  required 
all  visibility  predictions  be  doubled  before  applying  them  to 
Riley's  equation.  This  practice  was  selected  because  at  worst, 
it  errs  on  the  side  of  safety  in  any  thermal  calculation.  From 
an  aspect  of  troop  safety,  this  method  would  provide  for  more 
thermal  energy  at  the  receiver  location  due  to  an  inc  ease  in 
visual  range,  thereby  providing  less  attenuation  by  the 
atmosphere.  However,  it  does  introduce  some  uncertainty  into  the 
equation. 
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SECTION  3 


KSC  THERMAL  PREDICTION  METHOD 


The  KSC  Thermal  Prediction  Method  to  generate  the  probabil¬ 
istic  exposure  predictions  is  summarized  in  Figure  9.  A  crucial 
element  is  the  prediction  of  the  nuclear  source  output  which  is 
accomplished  by  the  RECIPE  code.  The  thermal  radiation  from  a 
nuclear  burst  was  not  explicitly  discussed  in  the  UK  documents. 
The  RECIPE  code  was  not  modified  on  this  contract,  but  is 
discussed  for  completeness  in  the  following  subsections.  The 
main  calculational  chain  is  shown  in  the  center  of  the  chart  and 
included:  (1)  a  weather  analysis  for  the  particular  area  of 

interest;  (2)  identification  of  pertinent  meteorological 
parameters  to  describe  the  broad  range  of  weather  conditions  of 
interest  primarily  involving  the  visibility  and  albedo  surface 
specifications;  (3)  computation  of  transmission  prediction  curves 
for  these  parameters;  and  (4)  converting  the  results  to 
transmission  occurrence  probabilities.  These  are  then  combined 
with  the  RECIPE  source  to  obtain  the  probabilistic  exposure 
predictions.  An  important  ingredient  in  the  TAXV  predictions  is 
a  very  large  data  base  of  build  up  factors  generated  by  the  KSC 
TRAX  Monte  Carlo  Code.  This  data  base  is  generated  for  a  wide 
range  of  meteorological  conditions  and  can  be  used  for  many 
diverse  applications.  If  exposure  predictions  are  desired  for  a 
specific  set  of  weather  conditions  then  the  probabilistic 
routines  are  bypassed  as  shown.  The  irradiance  at  a  receiver  is 
given  by: 

H  =  —2—  JT(  A,  Hs,  Hr,  R)  *P  (  A  ,  t,  Hs)  *  Eft(A,  Hs,  Hr,  R) 

4  JT  Ra^ 


where 


(7) 
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and 


P 

T 

X 


is  the  spectral  power  discussed  in  Section  3.1.1 

is  the  transmission  factor  discussed  in  Section  3.4 

is  the  wavelength 

is  the  source  altitude 

is  the  receiver  altitude 

is  the  slant  range 

is  the  time 

is  the  albedo  surface  enhancement  factor 
is  a  geometry  factor 


The  exposure  is  obtained  by  integrating  this  expression  over 
time  in  the  TAXV  code. 


The  transmission  is  given  by  the  product  of  the  direct 

transmission  (T  )  and  a  scattering  buildup  factor  (B  )  . 

D  S 

T(\,  Hs,  Hr,  R)  «  Td(X(Hs,  Hr,  R)  ‘Bs  (Ji,  Hs,  Hr(Ts)  (8) 

where  all  the  variables  have  been  defined  except  for  Tg  which  is 
the  scattering  optical  depth  found  by  integrating  the  scattering 
cross  section  along  the  sight  path  from  source  to  receiver. 

Implicit  in  these  expressions  is  a  model  atmosphere  containing  a 
particular  definition  of  atmospheric  parameters  as  a  function  of 
altitude  including  the  surface  visibility,  and  water  vapor,  and 
albedo  surface  definition. 


The  geometry  factor  and  the  power  are  computed  in  the  RECIPE 
module.  The  build-up  factor  and  the  albedo  enhancement  factor 
form  the  Monte  Carlo  data  base  and  are  computed  by  the  TRAX  Code. 
The  other  required  factors  and  computations  are  contained  in  the 
TAXV  modules.  In  the  following  sections  the  variation  of  the 
computational  scheme  are  described  in  detail. 
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3.1  RECIPE  SOURCE  MODULE. 


Only  a  minor  effort  on  this  contract  involved  the  RECIPE 
code  but  for  completeness  the  code  will  be  briefly  described 
since  it  forms  a  crucial  element  in  our  calculational  scheme  and 
will  be  used  in  future  thermal  predictions.  The  code  was 
originally  developed  by  SAI  (References  9  and  10)  over  a  period 
of  many  years.  Reference  11  utilizes  this  code  as  the  basis  for 
the  fireball  output  specifications.  The  code  has  undergone 
considerable  revision  and  simplification  at  KSC  (References  12 
and  13)  while  retaining  the  full  capability  of  predicting  the 
detailed  source  output  for  a  very  wide  range  of  burst  yields  and 
altitudes  including  surface  interaction  effects.  The  material 
given  here  is  extracted  primarily  from  the  KSC  portion  of 
Reference  13  and  is  included  here  because  the  balance  of 
Reference  13  will  be  of  no  interest  to  users  of  this  document. 

3.1.1  RECIPE  Code  Description. 

RECIPE  is  a  fast  running  Fortran  computer  code  that 
calculates  the  power  output  of  nucle  ir  fireballs  as  a  function  of 
time  and  wavelength.  The  original  version  contained  atmospheric 
transmission  and  instrumentation  routines  which  were  used  to 
analyze  the  nuclear  test  data.  The  unclassified  version  of 
RECIPE  was  used  for  these  studies.  A  classified  version  exists 
which  accurately  predicts  the  details  of  the  first  pulse 
signature  for  use  in  nuclear  detection  diagnostic  purposes.  The 
details  of  the  first  pulse  are  not  important  for  typical  damage 
response  analyses  since  essentially  all  of  the  energy  is 
contained  in  the  second  thermal  pulse.  For  detailed  analyses  of 
the  response  of  personnel  to  f lashblindnets  or  retinal  burns  and 
for  evaluation  of  goggles  or  photochromic  protective  devices,  a 
more  detailed  model  of  the  first  pulse  might  be  necessary. 
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The  original  version  developed  the  power  at  56  wavelengths 
between  200  and  12500  nm.  The  grid  of  56  wavelengths  was  chosen 
to  give  detailed  representation  of  the  atmospheric  transmission 
and  the  response  of  the  instrumentation  used  on  the  nuclear 
tests.  The  source  module  of  RECIPE  actually  uses  only  26 
wavelengths  to  describe  the  spectral  output  of  the  fireball.  To 
save  running  time  without  sacrificing  any  output  details,  the 
spectral  grid  of  the  retained  routines  was  reduced  to  the  basic 
26  groups. 


The  input  parameters  required  by  the  source  routines  of 
RECIPE  and  their  intervals  of  validity  are: 


Yield 
Altitude 
Bomb  Mass 


1  . E-4  to  3.E4  (kt) 

<30  (km) 

<2.E6  (lb/kt) 


The  power  output  is  given  as  a  function  of: 


Time 

Wavelength 
Elevation  Angle 


l.E-6  to  30*T2MAX  (sec) 

200  to  12500  (nm) 

-90  to  90  (degrees) 


The  evaluation  angle  dependence  arises  only  for  those  cases 
involving  surface  interactions  where  the  fireball  is 
rionspherical.  T2MAX  is  used  throughout  this  section  to  refer  to 
the  time  of  the  peak  of  the  main  or  second  thermal  pulse. 

The  RECIPE  source  module  contains  a  free  air  burst  model,  a 
surface  burst  model,  an  interacting  burst  model,  and  a  high 
altitude  burst  model.  These  models  parallel  the  historical 
development  of  the  code.  A  short  description  of  these  models 
will  be  given  including  the  routines  that  are  used  in  each  model. 
Additional  details  are  contained  in  References  9  and  10. 
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3. 1.1.1  Free  Air  Burst  Model. 


The  basic  data  used  in  the  free  air  burst  model  are  based 
upon  a  set  of  four  detailed  theoretical  radiation-hydrodynamic 
code  calculations  run  with  the  RHGEN  code  at  yields  of  4,  14, 

110,  and  6600  kt.  The  first  three  correspond  to  three  U.S. 
nuclear  atmospheric  bursts  which  have  been  measured  and  analyzed 
in  great  detail.  The  RHGEN  code  contains  the  latest  air  opacity 
data,  a  realistic  debris  mixing  model  and  an  improved  technique 
for  shock  front  location  which  is  necessary  to  provide  good 
definition  of  the  early  time  pulse  characteristics.  The  debris 
mixing  model  was  found  to  be  necessary  to  provide  good 
comparisons  with  measured  second  maxima  times  and  powers. 

The  mixing  of  the  bomb  debris  in  the  fireball  results  in 
changes  in  the  radiative  properties  of  the  fireball,  but  affects 
the  thermodynamic  properties  very  little.  The  increased  opacity 
due  to  the  debris  contribution  seems  to  slow  the  radial  heat  flow 
near  the  time  of  minimum,  increases  the  pressure  in  the  fireball, 
causing  a  continued  expansion  with  a  resulting  delay  in  the  time 
of  second  maximum.  The  increased  opacity  causes  an  increase  in 
the  effective  fireball  radius,  a  decreast  in  the  effective 
temperature,  and  a  lower  power. 

Two  different  procedures  are  used  in  the  RECIPE  model  to 
represent  the  RHGEN  results.  Prior  to  the  time  of  breakaway, 
which  corresponds  to  the  time  at  which  the  shock  becomes  trans¬ 
parent,  the  spectral  power  of  the  fireball  is  calculated  by 
combining  the  geometrical  properties  of  the  shock  with  the 
spectral  brightness  of  the  shock.  The  time  of  breakaway  is  about 
one-third  that  of  minimum  so  that  a  very  small  amount  of  energy 
is  released  prior  to  this  time.  The  details  of  the  handling  of 
this  portion  of  the  output  are  of  little  significance  in  response 
calculations  especially  at  low  altitudes. 
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After  the  time  of  breakaway  the  radius  and  brightness 
temperature  are  different  for  each  wavelength.  The  RECIPE  code 
uses  a  curve-fitting  procedure  to  reproduce  the  spectral  power 
curves  of  RHGEN.  The  eight  parameters  needed  are  the  time  of 
breakaway,  time  of  minimum,  time  of  second  maximum,  ten  times  the 
time  of  second  maximum  and  the  spectral  powers  at  these  times. 
These  parameters  are  used  in  the  code  to  generate  the  spectral 
power  at  an  arbitrary  time. 

The  routines  used  in  the  free  air  burst  model  are: 

SPFLUX:  principal  routine  which  calculates  the  spectral 

power  at  the  fireball  surface  at  a  given  time 
utilizing  the  following  subroutines. 

FRONT:  provides  shock  front  properties  as  a  function  of 

yield,  burst  altitude,  and  bomb  mass. 

RHSF:  provides  the  spectral  brightness  of  the  shock 

front  as  a  function  of  the  shock  velocity  and 
altitude. 

PARLOW:  provides  parameters  for  low  altitude  bursts  for 

utilization  by  the  formulas  in  SPFLUX. 

For  many  USANCA  applications  the  burst  is  low  enough  that 
significant  interaction  of  the  fireball  with  the  ground  surface 
occurs  resulting  in  a  significant  change  in  the  radiating 
characteristics  of  the  fireball.  These  effects  have  been 
incorporated  into  surface  interaction  modules  which  have  been 
thoroughly  check  with  available  atmospheric  test  data  (Reference 
10)  . 
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3. 1.1. 2  Surface  Burst  Model. 


The  surface  burst  model  is  based  to  a  large  extent  on 
analysis  of  nuclear  test  results  because  of  the  difficulty  in 
representing  the  experimental  results  entirely  from  theoretical 
results  (Reference  10) .  For  times  through  minimum  the  fireball 
development  is  the  same  as  that  for  a  free  air  burst  but  with  a 
different  effective  yield.  The  parameters  needed  to  obtain  the 
spectral  power  at  the  later  times  were  determined  from  experi¬ 
mental  data. 


3. 1.1. 3  Interacting  Burst  Model. 

Photographic  data  were  used  extensively  to  model  the 
fireball  development  that  occurs  when  surface  interactions  occur. 
The  fireball  at  early  times  does  not  touch  the  ground  but  during 
the  growth  phase  does  intersect  the  ground.  Two  different 
regimes  are  suggested:  The  interaction  occurs  early  enough  that 
a  transition  to  a  surface  burst  occurs,  and  an  interaction  occurs 
late  enough  that  the  fireball  shape  is  perturbed  but  remains 
above  the  ground. 

3. 1.1. 4  Surface  Burst  Transition. 

If  the  fireball  contacts  the  surface  prior  to  shock 
formation,  the  low  altitude  air  burst  is  transformed  into  a 
surface  burst.  The  spectral  power  computed  in  the  code  always 
refers  to  what  radiated  from  a  sphere.  For  the  surface  burst  an 
effective  yield  factor  of  1.96  is  used  and  a  shape  factor  (SAF) 
is  applied  to  represent  the  dependence  of  the  fluence  on  the 
elevation  angle  to  the  receiver. 

Figure  10  shows  how  this  is  modeled  as  a  function  of  time. 

In  the  top  illustration  for  the  air  burst  at  altitude  H(O)  the 
ratio  (WFAC)  of  the  effective  yield  (W)  to  the  device  yield  (WKT) 
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Figure  10.  Surface  Burst  Model. 


is  equal  to  one  at  the  early  times.  WFAC  is  computed  in  the 
following  manner  for  the  transition  positions  as  represented  by 
the  middle  illustration.  If  the  surface  were  not  there,  the 
fireball  wruld  develop  as  shown  by  the  dashed  lines.  Represent 
the  volume  above  the  surface  by  V  and  the  volume  below  the 
surface  by  V^.  If  the  assumption  is  made  that  a  fraction  R  of 
the  volume  below  the  surface  is  returned  to  the  above  ground 
perturbed  fireball  V',  that  the  centers  of  the  volumes  are  the 

cl 

same,  and  that  the  effective  yield  is  proportional  to  the  total 
volume,  then 


Va  ■  V  +  R*V. 

a  d 

(9) 

WFAC  =  v'  /V  =  v'/V  «  1  -  R  +  R/f 

cl  cl 

(10) 

where 

f  *  V/(V  +  v.  ) 
a  a  d 

(ID 

is  the  fractional  volume  above  the  surface. 

The  model  reduces  the  effective  burst  altitude  linearly  in 
time  so  that  the  burst  altitude  is  zero  at  twice  the  time  of 
first  contact  causing  the  sphere  to  transform  smoothly  into  a 
hemisphere.  The  assumed  value  of  R  is  .96  and,  therefore,  the 
limiting  value  for  WFAC  is  1.96. 

3. 1.1.5  Interacting  Air  Burst. 

For  this  case,  contact  with  the  surface  is  made  after  shock 
wave  formation,  and  the  assumption  is  that  negligible  energy  is 
transmitted  to  the  surface.  Figure  11  illustrates  the  model  used 
when  the  contact  occurs  prior  to  the  time  of  breakaway  as 
indicated  by  photographic  analysis.  The  truncated  sphere  is 
defined  in  terms  of  the  parameter  b  given  by  the  expression: 
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Figure  11.  Interfacing  Air  Burst  Model 


b ( t)  =1.24  (1  -  exp(-g) ) 


(radians) 


(12) 


where 

g  =  0.4  (t/te)**2.4  (13) 

t  =  time 

te  •  time  of  contact 

The  maximum  value  of  b  allowed  is  71  degrees.  The  volume  of  the 
truncated  sphere  is  equated  to  the  underformed  sphere. 

The  algorithm  for  the  case  where  the  shock  contacts  the 
ground  after  breakaway  is  similar  to  the  above  except  that  the 
fireball  radius  and  shock  front  radius  are  taken  into  account 
independently.  The  altitude  of  the  reflected  shock  is  computed, 
and  the  time  of  contact  for  use  in  the  above  formulas  is  taken  as 
the  time  at  which  the  reflected  shock  reaches  the  bottom  of  the 
fireball.  A  detailed  discussion  of  these  factors  is  contained  in 
Reference  10.  Since  the  fireballs  are  not  in  general  spherical 
the  shape  factor  depends  upon  the  elevation  angle  between  the 
fireball  symmetry  axis  and  the  receiver. 

The  routines  used  to  describe  the  interacting  bursts  are: 

WEFF:  main  routine  in  the  interacting  burst  model.  Sets 

interaction  times.  Provides  transition  in 
properties. 

SHAPE:  provides  shock  radius,  glowing  radius,  and  shock 

temperature. 

DARKR:  accounts  for  dirt  plume  for  surface  sources 
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3. 1.1. 6  High  Altitude  Model. 


The  low  altitude  models  were  extended  by  altitude  scaling 
techniques  valid  to  altitudes  of  30  km.  The  data  for  the 
altitude  scaling  of  the  shock  properties  needed  for  times  less 
than  breakaway  were  obtained  from  FIREBALL  runs. 

After  breakaway  the  same  type  of  parametric  fits  used  in  the 
free  air  model  are  used  with  altitude  scaling  using  the  ratio 
D  *  1.22E-3/p  where  p  is  the  ambient  density  at  a  selected 
altitude.  The  power  level  for  each  wavelength  is  given  by 

Pw  =  PQ  exp  (-r)  (14) 

where  Pw  is  the  power  ltvel  at  time  of  minimum  with  NC^  absorp¬ 
tion,  PQ  is  the  same  quantity  without  N02  absorption,  and  T 
represents  the  effective  optical  thickness  of  the  N02  for  a  given 
yield. 

Values  of  P  .  and  P^  are  available  at  sea  level  (SL)  and  thus 
w  o 

define  T.  At  a  selected  altitude  the  T  are  given  by 

t(D)  =  r(SL)  /D**l . 26  (15) 

at  each  wavelength. 

To  summarize,  the  set  of  RECIPE  routines  when  accessed  at  a 
certain  time  with  a  particular  yield,  altitude  and  bomb  mass 
return  the  fireball  output  power  for  each  of  the  26  wavelengths 
at  the  time  desired.  The  basic  calculations  give  the  output  for 
a  spherical  fireball,  and  a  shape  factor  is  applied  to  obtain  the 
power  in  a  particular  receiver  direction  when  that  information  is 
supplied.  The  shape  of  the  fireball  is  computed  at  each  time, 
an4,  a  <y  modifying  factors  are  computed  at  each  time  step.  For 
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incorporation  into  TAXV  the  possibility  of  a  rising  fireball  was 
included  in  the  routines,  so  that  the  shape  factors  computed  by 
the  code  includes  this  effect  which  was  neglected  in  the  original 
RECIPE  code. 

Figure  12  shows  the  thermal  fraction  as  a  function  of  yield 
with  the  burst  altitude  as  a  parameter.  This  represents  the 
fraction  of  the  total  device  yield  that  has  been  radiated  from 
the  fireball  up  to  10*T2MAX. 

A  general  downward  trend  is  noted  with  increasing  yield. 

The  basis  for  these  data  are  from  rad-hydro  runs  at  4,  14,  110, 
and  6600  kt  yields  with  the  yields'  dependence  being  determined 
by  the  interpolation  scheme  used  in  the  code. 

Figure  13  shows  the  same  information  plotted  as  function  of 
burst  altitude  with  the  yield  as  a  parameter.  These  curves  do 
not  include  the  reduction  due  to  surface  interaction  effects.  In 
Table  6  the  thermal  fractions  for  surface  and  free  air  bursts  are 
shown  along  with  the  transition  height.  Note  that  fr>r  the  higher 
yields  the  reduction  in  the  partition  for  the  surface  burst  is 
about  a  factor  of  two.  For  the  lower  .yields  the  reduction  in 
partition  is  much  larger. 
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080 


Figure  12.  Thermal  Yield  Fraction  as  a  Function  of  Burst 
Altitude  and  Yield  (Yield  Contours) 

(Source:  Hillendahl,  1980;  EM-1  Date:  February  1982) 
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30 


Thermal  Partition  (Etot/W) 

Figure  13.  Thermal  Yield  Fraction  as  a  Function  of  Burst 
Altitude  and  Yield  (Altitude  Contours) 

(Source:  Hillendahl,  1980:  EM-1  Date:  February  1982) 
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Table  6.  Thermal  Partition  for  Near-Surface  Bursts. 


Surface  Burst 

Non-Surface  Burst 

Transition 

Yield 

Partition 

Partition 

Height 

(KT) 

(Fraction) 

(Fraction) 

(Meters) 

1 

0.045 

0.350 

4 

10 

0.066 

0.341 

8.6 

100 

0.13 

0.330 

18.5 

1000 

0.16 

0.291 

40 

10000 

0.17 

0.254 

86 

The  difference  noted  for  the  surface  and  non-surface  bursts 
have  major  implications  in  predictions  of  the  exposure  . rom  low 
altitude  bursts  especially  for  receivers  near  the  ground  surface. 
For  lower  yield  tactical  devices  the  thermal  output  of  a  surface 
burst  is  ]  /3  that  of  a  low  altitude  free  air  burst  for  a  100  kt 
yield.  For  a  1  kt  yield  the  ratio  is  about  1/8.  A  shape  factor 
must  also  be  considered  since  the  fireballs  for  the  interacting 
bursts  are  strongly  perturbed  and  are  hemispherical  in  shape. 
Examples  of  the  magnitude  of  these  effects  on  the  exposure 
predictions  will  be  considered  in  the  following  sections. 

3.1.2  Time  Dependent  Power. 

Previous  USANCA  thermal  prediction  methods  have  ignored  time 
dependent  effects.  For  many  applications  and  certainly  for  any 
detailed  response  calculations  the  time  dependence  of  the  thermal 
environment  is  very  important.  For  eye  damage  effects  and 
especially  for  evaluation  of  eye  protective  devices,  the  time 
dependence  of  the  radiant  level  is  very  important.  The  RECIPE 
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code  gives  the  time  dependent  fireball  power  output  in  detail 
including  a  general  representation  of  the  first  pulse  as 
discussed  previously,  a  classified  version  is  available  giving 
the  first  pulse  in  detail. 

The  basic  output  of  the  SPFLUX  routine  is  the  spectral  power 
FOLZ(H  )  (w/eV)  as  a  function  of  photon  energy  at  the  time  of 
interest.  The  shape  factor  (SAF)  is  then  applied  to  obtain  the 
power  emitted  in  a  particular  receiver  direction 

FOLZH  (h v)  *  FOLZ(hp)  *  SAF.  (16) 

Two  alternate  time  dependent  modes  are  available.  In  one  mode  a 
time  mesh  is  defined  and  the  above  spectral  power  is  stored  for 
each  of  the  time  steps  for  use  in  the  predictive  routines  of 
TAXV.  More  detail  about  this  mode  will  be  given  later.  In  the 
other  mode,  the  calculation  above  is  completed  for  each  single 
time  step  defined  by  the  TAXV  control  routine  of  interest. 

The  power  is  then  obtained  as  a  function  of  time  by 
summation  over  the  energy  intervals 

P(t)  »  2F0LZH(hj/)  *  Ahv  .  (17) 

hi> 

t 

As  will  be  discussed  in  the  following  section  the  spectral  power 
is  actually  converted  to  a  wavelength  dependence  prior  to 
regrouping  into  the  wavelength  mesh  used  in  the  transmission 
routines  and  further  processed  into  the  standard  TAXV  wavelength 
grid.  In  Figure  14  are  shown  the  powertime  curves  for  100  kt 
burst  for  various  source  altitudes.  As  discussed  previously  the 
unclassified  version  of  RECIPE  has  been  used  in  this  development 
so  the  first  pulse  represents  only  the  continuum  contribution 
from  the  heated  shock.  The  details  of  the  NC^  absorption  and  its 
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Figure  14.  Effect  of  Altitude  on  Total  Thermal  Power, 
100-Kiloton  Burst. 
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effect  on  the  early  time  signature  are  not  included.  This  is  not 
important  for  material  response  considerations  since  only  a  very 
small  fraction  of  the  energy  is  contained  in  the  early  pulse. 

These  curves  show  the  expected  trends  with  altitude.  The 
power  level  of  the  first  pulse  increases  with  increasing 
altitude,  the  minimum  is  shallower  as  the  altitude  is  increased, 
and  the  second  pulse  becomes  narrower  and  higher  as  the  altitude 
increases.  Note,  the  relatively  large  difference  between  the 
free  air  burst  at  1  km  and  the  surface  burst.  The  surface  burst 
has  a  much  smaller  second  thermal  power  maximum  and  a  somewhat 
longer  time  to  second  maximum,  T2MAX.  The  first  pulse  for  the 
two  cases  is  seen  to  be  very  similar.  These  curves  represent  the 
power  leaving  the  fireball  surface.  In  order  to  compute  the 
irradiance  at  a  receiver  it  is  necessary  to  weight  the  FOLZS  by 
the  atmospheric  transmission  and  the  spectral  dependent  response 
function  before  summing  over  the  wavelength  bands.  This  will  be 
discussed  in  later  sections. 

This  option  of  computing  the  response  involves  accessing  the 
RECIPE  routines  at  each  of  the  times  required  in  the  TAXV  control 
routines  which  can  become  a  relatively  time  consuming  process. 
Another  option  mentioned  earlier  involves  generating  the  matrix 
FOLZH  (hi;,  t)  for  the  26  energy  mesh  points  and  a  predetermined 
time  grid.  The  total  exposure  is  then  found  by  integrating  over 
the  time  grid. 

The  time  mesh  is  determined  in  the  following  manner.  A 
total  of  95  time  factors  are  defined  with  60  factors  increasing 
in  increments  from  .05  to  3.,  then  with  35  factors  increasing  in 
increments  of  .2  from  3.2  to  10.  The  actual  times  are  then  found 
by  multiplying  the  time  factors  by  the  time  of  second  maximum. 

In  this  manner  fine  time  steps  are  defined  from  zero  to  3  times 
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the  time  of  second  maximum  where  the  power  is  a  rapidly  varying 
function  of  time.  The  log-log  scales  in  Figure  14  do  not  give  a 
good  feel  for  the  actual  time  dependence  of  the  power  curves  for 
the  second  pulse.  In  Figure  15  the  power-time  curves  in  Figure 
14  are  replotted  on  normalized  linear  scales.  The  abscissa  is 
the  ratio  of  the  time  to  T2MAX.  The  second  thermal  pulse  when 
plotted  in  this  normalized  manner  is  seen  to  display  a  relatively 
small  altitude  dependence.  The  increase  in  the  minimum  power  at 
higher  altitudes  is  shown,  and  the  relative  unimportance  of  the 
first  pulse  in  terms  of  total  power  is  shown. 

In  Figure  16  the  irradiance  is  shown  for  bursts  at  an 
altitude  of  1  kft.  The  receiver  is  at  an  altitude  of  1  kft  and 
at  a  range  of  31  kft  for  the  lMt  and  at  21  kft  for  the  100  kt. 

The  total  exposure  for  both  cases  is  20  cal/cma.  The  data  points 
are  the  times  at  which  TAXV  computes  the  power  time  mesh  and  are 
seen  to  represent  the  shape  of  the  curves  in  fine  detail. 

The  code  does  not  contain  an  explicit  formula  for  the  total 
power  integrated  over  the  spectrum  at  T2MAX  nor  for  T2MAX  itself. 
Instead  curve  fits  are  used  in  the  code  for  each  of  the  26  energy 
grid  points.  The  expression  is  of  the  form: 

T2MAX  (h  v)  =  3.682E-2  *  (W**C(hp))  *  [p/pQ)**.  315  (18) 

where  the  parameter  C  is  a  function  of  hv,  and  the  altitude 
dependence  is  given  by  the  density  ratio  expression.  The  above 
formula  is  used  if  bomb  mass  is  less  than  2.5E3  lb/kt.  For 
heavier  bombs,  a  factor  involving  bomb  mass  is  included  which 
increases  T2MAX .  In  the  routines  developed  for  this  program 
T2MAX  for  the  total  integrated  power  is  represented  by  T2MAX  for 
550  nm  since  the  spectrum  tends  to  be  peaked  in  the  visible 
portion  of  the  spectrum. 
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Figure  16.  Irradiance  Predicted  with  the  Standard  Time  Factor 
Mesh  for  Bursts  at  1  KFT  Altitude  Coaltitude  at  the 
20  Cal/CM2  Exposure  Level. 
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3.1.3  Spectral  Dependence. 


Previous  USANCA  prediction  techniques  have  not  explicitly 
addressed  the  spectral  dependence  of  the  thermal  energy.  The 
spectral  distribution  of  the  power  from  the  fireball  is  of 
importance  primarily  for  determining  the  atmospheric  transmission 
from  the  fireball  to  the  receiver  and  also  to  a  lesser  extent  in 
determining  the  absorptivity  of  the  material. 

As  was  discussed  in  the  last  section,  the  spectral 
dependence  of  the  power  in  RECIPE  is  described  by  26  energy 
groups  which  span  the  wavelength  range  from  200  to  12500  nm.  All 
of  the  curve  fitting  in  the  code  is  done  as  a  function  of  these 
energy  groups  and  the  basic  power  matrices  are  developed  with  the 
units  of  watts  per  eV.  These  spectral  power  matrices  are  used  to 
compute  the  energy  in  calories  per  second  in  each  wavelength  band 
before  being  combined  with  the  atmospheric  transmission  factor 
and  the  spectral  response  functions  in  order  to  calculate  the 
total  energy  deposited  in  the  material  of  interest. 

In  Figure  17  the  spectral  distribution  is  shown  for  100  kt 
burst  at  several  burst  altitudes.  These  curves  are  obtained  by 
integrating  the  power  matrices  over  time  and  represent  the 
distribution  for  the  total  energy  radiated  from  the  fireball. 

Each  curve  has  been  normalized  to  unity.  The  free  air  bursts  all 
have  essentially  the  same  distribution  except  the  bursts  at 
higher  altitudes  tend  to  be  more  sharply  peaked  and  contain  more 
energy  in  the  UV  portion  of  the  spectrum.  The  surface  burst  is 
seen  to  be  definitely  shifted  to  longer  wavelengths  indicating 
the  effects  of  the  ground  surface  on  reducing  the  radiating 
temperature  of  the  fireball  and  the  increased  absorption  from  the 
entrained  material.  Most  of  the  energy  is  concentrated  in  the 
visible  portion  of  the  spectrum  with  a  small  fraction  having  a 
wavelength  greater  than  one  micron. 
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Source  Distribution 


Figure  17.  Effect  of  Altitud 
100  Kilotons. 
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In  Figure  18  the  source  distribution  is  displayed  in  a 
different  manner  for  100  kt  bursts  at  a  1  km  altitude  and  on  the 
surface.  The  fraction  of  the  energy  in  the  TRAX  (KSC  Monte  Carlo 
transmission  code)  wavelength  intervals  is  shown  for  the  two 
different  altitudes.  The  surface  burst  definitely  contains  a 
larger  fraction  of  the  energy  at  the  longer  wavelengths.  As 
shown,  the  QRT  data  base  wavelengths  grid  is  not  sufficient  to 
adequately  describe  the  spectral  output.  This  will  be  discussed 
in  Section  5. 

In  Figure  19  the  spectral  distribution  for  times  of  . 5T2MAX, 
T2MAX  and  10T2MAX  as  well  as  the  integration  over  time  is  shown 
for  the  100  kt  example  considered  in  the  previous  section.  There 
is  not  a  strong  time  dependence  in  the  spectrum,  but  the  full 
detail  is  available  in  the  RECIPE  code  and  is  utilized  in  the 
updated  routines.  Similar  spectral  distributions  are  shown  in 
Figure  20  for  the  1  Mt  case.  These  types  of  variations  will  not 
change  the  time  dependent  energy  absorptions  by  a  large  amount 
unless  the  absorptivity  shows  a  very  large  wavelength  dependence. 
Analyses  involving  windshields  and  specialized  materials  such  as 
photochromic  panels  will  require  this  type  of  detail. 

3.1.4  Fireball  Radius  and  Altitude  Model. 

Routines  have  been  added  to  TAXV  which  were  derived  from 
RANCIV  routines  (Reference  15)  which  give  the  fireball  center  as 
a  function  of  altitude.  The  original  RANCIV  routines  give  both 
the  altitude  and  dimensions  of  the  fireball  as  a  function  of  time 
for  particular  burst  altitudes  and  yields.  The  progression  of 
the  fireball  from  a  sphere  to  an  oblate  spheroid  to  a  torus  is 
represented.  The  later  time  development  into  a  torus  as  modeled 
by  RANCIV  is  not  of  importance  for  thermal  applications  because 
the  radiant  output  is  much  less  than  at  earlier  times.  The 
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Figure  18.  Source  Distribution  for  TRAX  and  QRT  Codes. 
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Figure  19.  Spectral  Distribution  Received  a  Distance  of  21 
KFT  Coaltitude  From  a  100  KT  Burst. 
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Figure  20.  Spectral  Distribution  Received  a  Distance  of 
31  KFT  Coaltitude  From  a  1  MT  Burst. 
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combination  of  the  RECIPE  dimension  modeling  including  the  early 
time  distortions  due  to  surface  effects  and  the  simplified  RANCIV 
fireball  altitude  modeling  results  in  simple  fast  running 
algorithms  adequate  for  thermal  response  calculations. 

Especially  considering  the  fact  that  for  USANCA  applications,  the 
receiver  point  is  seldom  close  to  the  fireball.  The  NORSE/ROSCOE 
routines  under  current  development  were  not  used  because  versions 
are  not  available  for  distribution  and  because  the  complexity  of 
these  routines  including  multiburst  effects  is  not  warranted. 

3.2  KSC  MONTE  CARLO  DATA  BASE. 


The  TRAX  code  was  developed  at  KSC  to  generate  a  Monte  Carlo 
data  base  of  thermal  transport  results  in  the  atmosphere,  speci¬ 
fically  aimed  at  nuclear  weapons  effects  applications.  The  TAXV 
code  was  developed  to  utilize  the  TRAX  data  base  in  the  most 
efficient  and  accurate  manner  and  to  easily  couple  to  the  RECIPE 
source  output  model  for  thermal  exposure  calculations.  These 
codes  have  been  exercised  and  improved  by  KSC  as  a  result  of 
several  DNA,  BMD,  and  USANCA  contracts. 

TRAX  is  a  very  general  Monte  Carlo  transport  code  that  was 
developed  specifically  for  generating  thermal  exposure  predic¬ 
tions  from  nuclear  weapons  in  the  atmosphere.  The  code  was 
originally  developed  during  the  preparation  of  Reference  12  which 
contains  a  comprehensive  discussion  of  all  thermal  data  collected 
on  the  U.S.  nuclear  test  series,  a  theoretical  description  of 
thermal  output  from  fireballs,  a  detailed  discussion  of 
atmospheric  transmission  concentrating  on  applications  to  thermal 
exposure  predictions,  and  an  empirical  derivation  of  the  thermal 
partition  as  a  function  of  yield  and  altitude.  The  TRAX  code 
utilized  previous  work  done  by  R1 3  in  atmospheric  transmission, 
but  incorporated  many  features  that  corrected  shortcomings  in 
their  approach  which  are  still  contained  in  the  QRT  code  (Refer¬ 
ence  14)  .  The  TRAX  code  was  developed  with  very  efficient  path 
selection  techniques,  efficient  broad  band  absorption  techniques, 
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and  a  2-D  scoring  technique.  The  scoring  technique  uses  a  path 
extension  approach  which  increases  scoring  efficiency  with 
minimal  computing  time  increase.  The  RRA  codes  rely  heavily  on 
theoretical  narrow  band  IR  absorption  techniques  while  the  TRAX 
code  uses  broad  band  experimental  results  to  most  efficiently 
handle  the  energy  in  the  IR  portion  of  the  nuclear  output 
spectrum. 


Eighteen  wavelength  bands  are  currently  used  in  the  calcu¬ 
lations  being  done  at  KSC.  These  bands  are  indicated  in 
Table  7  with  the  absorption  type  also  shown  as  originally  defined 
for  the  TRAX  code.  For  each  band  in  addition  to  the  absorption 
reactions  shown,  the  Rayleigh  scattering  from  the  molecular 
constituents  and  Mie  scattering  from  the  aerosol  particles 
suspended  in  the  atmosphere  are  considered.  The  cross  sections 
and  the  angular  phase  functions  for  both  scattering  functions  are 
calculated  at  the  midpoint  of  the  bands.  As  shown  in  section 
3.1.3  this  band  structure  provides  good  coverage  of  the  nuclear 
output  spectrum. 


The  transmission  (T^  )  along  a  path  is  computed  by  the 
expression: 


TX 


*  e 


'  V *  (X*hl}~  (X'Vl  X  R  X  TIR(X,hlfh2,R)  (19) 


hl  "  h2 


where  h^,  h2 


altitudes  at  the  extremes  of  the  path 
segment 


r^X,  h^)  *  total  extinction  optical  depth 
R  *  slant  path  length 
X  =  wavelength 


T  (X,h‘,h_,R)  =  IR  absorption  factor. 

IK  1  Z 
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If  h^  *  then  the  argument  of  the  exponential  is  replaced 

by  the  quantity  ae  (X#h^)R  where  0&  (X#^)  is  the  total 
extinction  cross  section  which  is  given  as  the  sum  of  the 
Rayleigh  and  Mie  scattering  cross  section  and  the  ozone  absorp¬ 
tion  cross  section. 

Table  7.  Trax  Band  Structure. 


BAND  WAVELENGTH 

<M> 

BAND  LIMITS 

/.LSORPTION 

TYPE 

0.25 

0.2  - 

0.3 

Ozone 

0.32 

0.3  - 

0.34 

Ozone 

0.36 

0.34  - 

0.4 

Ozone 

0.45 

0.4  - 

0.5 

Ozone 

0.55 

0.5  - 

0.6 

Ozone 

0.65 

0.6  - 

0.7 

Ozone 

3.80 

0.7  - 

0.8696 

Ozone 

0.94 

0.8696  - 

0.9901 

h2o 

1.03 

0.9901  - 

1.075 

None 

1.10 

1.075  - 

1.205 

h2o 

1.23 

1.205  - 

1.253 

None 

1.38 

1.253  - 

1.538 

h2o  + 

co2 

1.60 

1.538  - 

1.695 

co2 

1.87 

1.695  - 

2.083 

h2o  + 

C°2 

2.17 

2.083  - 

2.273 

None 

2.70 

2.273  - 

2.994 

h2o  + 

C02 

3.20 

2.994  - 

3.571 

h2o 

3.80 

3.571  - 

4.00 

None 

The  altitude  mesh  used  was  1  km  steps  up  to  50  km  but  this 
mesh  is  easily  changed  and  extended.  For  Rayleigh  and  Mie 
scattering  and  ozone  absorption.  The  cross  section  °R,  0Q, 

and  the  optical  depth  fR,  and  are  given  as  a 

function  of  altitude  and  wavelength.  The  total  extinction  cross 
section  0  and  the  total  optical  depth  T  e  is  also  given.  The 
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optical  depths  are  defined  from  the  top  of  the  atmosphere  down¬ 
ward  instead  of  the  other  direction  because  less  precision  is 
required  to  accurately  compute  differences  needed  in  the  calcula¬ 
tions.  Defining  the  optical  depths  from  the  bottom  may  result  in 
conditions  where  the  inherent  precision  of  a  computer  may  be 
stressed  in  calculating  the  difference  in  the  optical  depths 
between  two  altitudes. 

The  IR  transmission  factor  is  calculated  using  experimental¬ 
ly  derived  parameters  for  the  defined  bands  in  the  following 
manner.  The  basic  data  give  the  band  absorption  (A)  as  a 
fur.-tion  of  the  amount  of  absorber  in  the  path  segment  and  the 
average  pressure  along  the  path.  For  small  amounts  of  absorber  a 
weak  fit  is  used: 

A  =  CwV1  (20) 

and  for  larger  amounts  of  absorber  a  strong  fit  is  used: 

A  '  C  +  D  log  W  +  K  log  P  (21) 

where  W  is  the  integrated  absorber  and  P  is  the  pressure.  A 
transition  absorption  is  defined  between  the  two  fits.  For  C02 
absorption  then  the  transmission  for  the  band  is  given  by: 


1  -  A, 


C02 


CO. 


A  v 


(22) 


where  A  p  is  the  band  width  in  wave  numbers.  The  parameters 
used  are  given  in  Table  8. 
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Table  8.  Carbon  Dioxide  Absorption  Parameters* . 


Band  (fx) 

Weak 

C 

Fit 

K 

Transition 

A 

C 

Strong 

D 

Fit 

K 

1.38 

.058 

.41 

80 

- 

- 

- 

1.60 

.063 

.38 

80 

- 

- 

- 

1.87 

.492 

.39 

80 

-536 

138 

114 

2.70 

3.150 

.43 

50 

-137 

77 

68 

For  water  vapor  absorption,  it  was  found  that  at  large 
amounts  of  absorber,  the  laboratory  measurements  did  not  agree 
with  measurements  performed  in  the  atmosphere  for  these  bands.  A 
transition  absorber  was  defined  beyond  which  a  power  law  fit  to 
the  transmission  was  used: 

TH20  -  AW6  (23) 

In  Table  9  the  parameters  used  for  water  vapor  absorption  are 
listed.  For  bands  in  which  both  water  vapor  and  C02  absorption 
occur  the  expression: 

A  “  AH20  +  (1  “  2,36  AH20)  AC02  (24) 

was  used.  Using  these  broad  band  expressions  rather  than  narrow 
band  values  to  represent  a  broad  band  provides  a  more  accurate 
description  of  the  IR  absorption.  The  above  absorption  model  has 
been  compared  with  that  obtained  by  integrating  the  Air  Force 
Geophysical  Laboratory  (AFGL)  narrow  band  absorption  over  the 
broad  bands  given  above  and  a  good  agreement  was  found. 

The  TRAX  code  generates  histories  in  a  standard  manner  using 
the  scattering  cross  sections  to  define  path  lengths,  the 

*  From  HBW,  Reference  15. 
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Table  9.  Water  Vapor  Absorption  Parameters 
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scattering  phase  functions  to  define  particle  directions,  and  at 
each  collision  adjusts  the  weight  of  the  particle  to  allow  for 
the  possibility  of  being  absorbed  along  the  path.  Because  of  the 
nonlinear  characteristics  of  the  broad  band  absorption,  the 
cumulative  absorber  depth  traversed  is  used  to  compute  the 
particle  weight.  A  form  of  extended  path  length  scoring  is  used 
which  reduces  the  variance  of  the  Monte  Carlo  score.  This  is 
done  by  computing  the  probability  that  a  path  segment  could  have 
reached  any  scoring  surface  in  the  direction  of  propagation  if 
the  scattering  had  not  taken  place.  This  provides  a  safe  method 
of  increasing  the  efficiency  of  the  code  without  introducing  any 
bias  in  the  answer  as  other  variance  reduction  techniques  do 
unless  care  :ully  used. 

The  scoring  is  performed  in  a  defined  altitude  and  hori¬ 
zontal  range  mesh  and  with  sampling  angles  defined  relative  to 
the  source  to  receiver  vector.  The  scoring  gives  the  scattered 
flux  received  on  a  flat  plate  detector  at  the  location  and  with 
the  direction  defined.  The  scattered  flux  is  itself  never  used 
directly  for  exposure  predictions.  Instead,  the  scattered  flux 
is  used  to  determine  the  buildup  factor  which  is  a  ratio  of  the 
total  flux  received  at  the  receiver  point  (direct  flux  plus  the 
scattered  flux)  divided  by  the  direct  flux. 

The  buildup  factor  shows  a  much  more  predictable  variation 
with  geometry  parameters  and  wavelength  than  does  the  scattered 
flux.  In  the  IR  region  as  shown  in  Table  7,  the  defined  bands 
include  bands  with  strong  absorption  and  also  bands  where  no 
absorption  occurs  representing  windows  in  the  spectrum.  In 
Figure  18,  typical  results  are  shown  for  an  atmosphere  defined 
for  Nevada  atmospheric  conditions  for  several  wavelengths.  In 
this  case  the  total  transmission  is  plotted  as  a  function  of 
horizontal  range.  The  .32  and  the  .55  bands  do  not  include 
absorption  and  the  lower  transmission  for  .32  is  due  to  the  much 
larger  scattering  cross  sections  at  the  lower  wavelengths.  In 
the  IR  region,  note  that  .94  shows  a  relatively  strong  absorp- 
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tion,  1.23  shows  essentially  no  absorption,  and  1.87  shows  strong 
absorption.  The  1.23  band  has  no  IR  absorption  and  the  reduction 
in  the  transmission  is  due  only  to  the  small  scattering  that 
occurs  at  this  wavelength.  A  plot  of  the  transmission  or  the 
total  flux  which  is  an  equivalent  quantity  as  a  function  of  the 
wavelength  would  show  no  obvious  trend  with  wavelength  as  will  be 
shown  in  Section  5  in  discussing  the  QRT  approach. 

However,  when  the  buildup  factors  are  plotted  as  a  function 
of  wavelength  as  is  done  in  Figure  22,  a  very  obvious  trend  is 
observed.  There  is  seen  to  be  a  very  well  defined  trend  in  the 
buildup  factor  versus  scattering  optical  depth  as  one  expects 
from  elementary  transport  theory,  but  also  a  well  defined  varia¬ 
tion  with  wavelength  which  might  not  have  been  predicted.  The 
cross  section  decreases  with  increasing  wavelength  in  a  smooth 
manner  and  likewise  the  buildup  factor  also  decreases  with 
wavelength.  The  magnitude  of  the  total  flux  at  a  point  in  space 
depends  strongly  on  the  wavelength,  but  the  ratio  of  the  total 
flux  to  the  direct  flux  does  not.  This  forms  the  basis  for  all 
of  the  predictive  codes  used  at  KSC  and  allows  for  much  more 
accurate  scaling  of  the  Monte  Carlo  results  than  is  obtained  by 
using  the  scattered  flux  as  is  done  in  the  QRT  code. 

A  problem  with  the  current  KSC  data  base  is  that  the  data 
were  developed  prior  to  a  redefinition  of  model  atmospheric 
parameters  at  AFGL.  The  newer  models  have  been  incorporated  into 
newer  TAXV  routines  as  described  in  Section  3.4. 

Sample  buildup  factor  results  have  been  generated  with  the 
TRAX  results  to  define  the  applicability  of  the  current  results. 
The  buildup  factors  for  the  data  base  aerosol  distribution  and 
the  rural  and  maritime  aerosol  distributions  are  shown  in  Figure 
23  as  a  function  of  scattering  optical  depth  for  a  wavelength  of 
.55  microns  for  a  meteorological  range  of  23  km.  The  1  sigma 
bounds  for  the  data  base  points  are  shown.  Minor  differences  are 
noted  for  scattering  optical  depths  of  about  1.3  corresponding  to 
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TOTAL  TRANSMISSION 
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Figure  21.  Comparison  of  Tax  Results  with  TRAX  Data  Nevada 
Atmosphere  with  Ground  Level  at  1.28  KM  Source 
1  KM  Above  Ground.  Receiver  3  M  Above  Ground. 
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TOTAL  RADIATION/  DIRECT  RADIATION 


Figure  22.  Comparison  of  Buiid-up  Factors  vs  the  Scattering 
Optical  Depth  for  Various  Wavelengths.  Pacific 
Atmosphere  with  Both  Source  and  Sampling  at 
1  km  Altitude. 
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Figure  23.  Comparison  of  Build-up  Factors  for  Surface  Burst. 
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Scattering  Optical  Depth 


ranges  of  7.5  km.  Figure  24  and  25  show  corresponding  results 
for  meteorological  ranges  of  10  and  2  km.  The  corresponding 
ranges  9.5  and  1.5  km  for  10  and  2  km  meteorological  ranges 
respectively  show  which  data  base  results  are  acceptable. 

For  100  kt  yields  the  predictions  are  acceptable  down  to 
2  cal/cm2  for  meteorological  ranges  greater  than  10  km.  For 
lower  meteorological  ranges  the  predictions  are  valid  for 
exposures  greater  than  about  40  cal/cma.  For  1  Mt  yields  the 
exposure  predictions  curves  are  good  for  exposures  greater  than 
20  cal/cm2.  For  low  meteorological  ranges  the  prediction  are 
valid  only  for  exposures  over  about  400  cal/cm2. 

These  boundary  conditions  indicate  that  data  base  revisions 
are  necessary  to  adequately  define  the  probabilistic 
transmissions  for  USANCA  applications.  As  discussed  above  the 
transmissions  are  valid  only  for  large  meteorological  ranges. 

The  probabilistic  transmissions  are  accurate  only  for  the  higher 
probabilities . 

3.3  WEATHER  ANALYSIS  REQUIREMENTS. 

In  the  previous  section  the  UK  weather  analysis  was 
discussed.  The  mai^  factors  considered  in  their  analysis  was  the 
visibility  and  the  cloud  ceiling.  Other  factors  need  to  be  con¬ 
sidered  in  order  to  perform  accurate  thermal  predictions. 

A  set  of  weather  parameters  is  needed  for  developing  the 
basic  Monte  Carlo  data  base.  As  discussed  in  the  previous 
subsection  a  very  detailed  set  of  atmospheric  parameters  are 
needed  for  these  calculations.  These  include  the  aerosol  and 
Rayleigh  cross  section,  water  vapor,  carbon  dioxide,  density 
ozone,  temperature  and  pressure  as  a  function  of  altitude.  Also 
needed,  are  the  Rayleigh  and  aerosol  scattering  phase  function. 
The  latter  of  which  depends  upon  the  specific  aerosol  size  dis- 
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BUILDUP  FACTOR 


Figure  24.  Comparison  of  Build-up  Factors  for  Surface  Burst. 
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Figure  25. 


tribution  chosen.  The  ground  surface  and  cloud  cover  altitudes 
and  reflecting  characteristics  are  also  needed. 

The  basic  atmospheric  constituents  chosen  for  future  work 
correspond  to  the  mid-latitude  fall-winter  and  the  mid-latitude 
spring-summer  conditions  (Reference  17)  .  The  spring-summer  model 
has  a  higher  water  vapor  concentration  than  the  fall-winter  model 
and  represents  the  constituent  having  the  widest  variation. 

The  aerosol  model  chosen  (Reference  14)  corresponds  to  a 
rural  aerosol  distribution  with  a  ground  surface  meteorological 
range  of  23  km.  The  basic  parameter  is  the  aerosol  cross  section 
(km  *)  for  the  .55  micron  band  as  a  runction  of  altitude.  This 
parameter  is  the  same  for  both  model  atmospheres  at  the  lower 
altitudes.  The  code  permits  a  different  meteorological  range  at 
the  surface  to  be  chosen,  and  the  aerosol  cross  section  is 
adjusted  at  the  lower  altitudes  in  the  boundary  layer  to  obtain 
the  required  meteorological  range.  A  similar  procedure  is  used 
to  alter  the  water  vapor  distribution  at  the  lower  altitudes. 

In  Table  10  the  aerosol,  Rayleigh  and  ozone  parameters  are 
given  for  the  spring-summer  model  for  the  .55  micron  band.  For 
each  constituent  both  the  cross  section  o(km  )  and  the  optical 
depth  r  are  given.  The  total  extinction  cross  section  and 
optical  depth  are  obtained  by  summing  the  three  components. 
Corresponding  values  for  the  fall-winter  model  are  given  in 
Table  14. 

Figures  26,  27  and  28  show  the  Rayleigh,  aerosol  and  ozone 
cross  sections  as  a  function  of  altitude.  The  solid  line  refers 
to  the  fall-winter  atmosphere  and  the  dashed  line  refers  to  the 
spring-summer  atmosphere. 

These  parameters  are  stored  in  the  code,  and  whenever  values 
of  <7e  or  Tg  are  needed  to  compute  the  transmission,  interpola- 
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Figure  26.  Comparison  of  Rayleigh  Cross  Section  for  Two  Model 
Atmosphere . 


70 


ALTITUDE  (KM) 


Figure  27.  Comparison  of  Aerosol  Cross  Section  for  Two  Model 
Atmosphere . 
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Figure  28. 


Comparison  of  Ozone  Cross  Section 
Atmosphere . 
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Scattering  Parameters  for  the  Spring-Summer  Model  Atmosphere 
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Table  11.  Scattering  Parameters  for  the  Fall-Winter  Model  Atmosphere 
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Table  12.  Wavelength  Dependent  Cross  Section  Factors 


WAVELENGTH 

BAND 

RAYLEIGH 

AEROSOL 

OZONE 

.25 

21.84 

1.77 

1.82E3 

.32 

9.48 

1.58 

9.76 

.36 

5.76 

1.52 

1.95E-2 

.45 

2.28 

1.14 

3.81E-2 

.55 

1.0 

1.0 

1.0 

.65 

5.08E-1 

0.899 

6.74E-1 

.80 

2.19E-1 

0.804 

1.09E-1 

.94 

1.15E-1 

0.747 

0.0 

1.03 

7.91E-2 

0.722 

0.0 

1.10 

6.05E-2 

0.709 

0.0 

1.23 

3.86E-2 

0.684 

0.0 

1.38 

2.45E-2 

0.665 

0.0 

1.6 

1 . 36E-2 

0.630 

0.0 

1.87 

7.23E-3 

0.595 

0.0 

2.17 

3.99E-3 

0.557 

0.0 

2.7 

1.67E-3 

0.506 

0.0 

3.2 

8.10E-4 

0.465 

0.0 

3.8 

4.03E-4 

0.415 

0.0 
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tion  in  these  tables  is  performed.  The  tables  contain  values  for 
the  .55  p  band.  Values  for  other  bands  are  calculated  using  the 
ratios  (K)  in  Table  12  as  follows: 


<M  X  ,H)  ■  K  (  A  )*  <J( .  55 ,  H) .  (25) 

The  absorption  expressions  utilize  the  integrated  amount  of 
absorber  (w)  along  the  source  to  receiver  path.  If  Hg  ^  Hr  the 
amount  of  absorber  is 

w  *  |  w  (Hg)  -  w(Hr)l  *R  |Hs  -  Hrl  (26) 

where 


and  p  is  the  absorber  concentration  (either  C02  or  H20) .  If  Hg  ■ 
Hr  then 

w  •  p  *  R.  (28) 
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For  horizontal  paths  the  pressure  (P)  used  in  the  expression  is 
the  pressure  at  the  relevant  altitude.  If  Hg  #  Hr  then  an 
average  pressure  is  used  defined  by 

Pw (HS)  ‘  Pw(Hr) 

P  *  -  (29) 

w(H  )  -  w (H  ) 
s  r 

where 

50 

P„'Hs’  *  { 


The  CO j  density  in  these  expressions  is  expressed  in  units 
of  atmosphere-cm/km  and  the  H20  density  is  expressed  in  units  of 
precipitable-cm/km.  The  unit  atmosphere-cm  refers  to  the  depth 
of  absorber  resulting  from  the  absorber  traversing  along  a  path 
which  was  reduced  to  STP  conditions.  This  quantity  is  computed 
by  integrating  the  absorber  density  along  the  path  and  dividing 
by  the  STP  density.  The  carbon  dioxide  concentration  is  assumed 
constant  at  0.0314  percent  by  volume  for  all  altitudes  with  the 
number  density  as  a  function  of  altitude  determined  by  the  model 
atmosphere  used.  The  units  pr-cm  are  analogous  with  the  defini¬ 
tion  of  atmosphere-cm  where  it  is  assumed  that  the  water  vapor  is 
converted  to  the  liquid  state  of  STP.  Since  the  density  at  STP 
is  essentially  unity,  numerically  is  just  the  integral  of 

the  water  vapor  is  g/cma  along  the  path. 

0  r  w  W 

The  quantities  1^0,  C02,  H20  and  CC^  as  a  function  of 

altitude  as  well  as  the  quantities  weighted  by  the  pressure  are 
needed  to  compute  the  IR  absorption.  In  Table  13  the  quantities 

/i  W 

^C02,  C02,  and  C02  are  given  as  a  function  of  altitude 


PpdH. 


(30) 
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for  the  spring-summer  atmosphere  model.  Corresponding  values  are 
given  in  Table  14  for  the  fall  winter  model  atmosphere.  Likewise, 
in  Tables  15  and  16  the  corresponding  values  are  given  for  the 
water  vapor  concentration. 

In  Figure  29  the  carbon  dioxide  concentrations  for  the  two 
model  atmospheres  are  compared  with  the  solid  line  representing 
the  spring-summer  and  the  dashed  line  representing  the  fall- 
winter  model.  In  Figure  30  the  water  vapor  density  for  the  two 
models  is  compared.  The  spring-summer  model  is  seen  to  be 
significantly  wetter  than  the  fall-winter  at  altitudes  below 
about  10  km.  The  units  are  pr-cm/km  and  the  values  must  be 
multiplied  by  10  to  obtain  g/m3. 

The  previously  described  atmospheric  parameters  are  required 
for  the  basic  Monte  Carlo  data  base  requirements  and  are  also 
required  as  a  basis  for  the  prediction  code  TAXV  to  be  described 
in  the  next  section.  In  a  specific  transmission  prediction  run 
the  weather  parameters  for  that  case  need  to  be  specified.  The 
primary  parameters  of  interest  are:  the  ground  and  cloud  ceiling 
altitudes  and  the  percentage  cloud  cover,  the  visibility  and  the 
water  vapor  density  along  the  ground  surface,  the  model  atmo¬ 
sphere  type,  and  the  ground  albedo. 

The  ground  surface  visibility  and  water  vapor  density  are 
used  to  define  atmospheric  properties  in  the  boundary  layer  near 
the  ground.  The  cloud  and  ground  parameters  are  used  to  define 
albedo  surface  enhancement  factors  in  the  calculations. 

In  the  UK  weather  analysis  cloud  ceiling  above  3  km  was 
equated  with  the  absence  of  clouds.  Actually,  the  transmission 
factor  for  a  cloud  ceiling  at  6  km,  for  instance,  is  much  larger 
than  that  obtained  with  no  cloud  layer.  The  95%  curve  obtained 
with  clouds  is  much  higher  than  the  curve  without  clouds.  The 
effect  this  assumption  had  on  thermal  transmissivity  curves  can 
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Table  13.  Carbon  Dioxide  Absorption  Parameters  for  the  Spring-Summer 
Model  Atmosphere. 
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Table  14.  Carbon  Dioxide  Absorption  Parameters  for  the  Fall-Winter 
Model  Atmosphere. 
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L5.  Water  Vapor  Absorption  Parameters  For  The  Spring-Summer  Model  Atmosphere 

WATER  VAPOR  _  PRESSURE  *  WATER  VAPOR 

DfcXSITY  DEPTH  DEXSITY  DEPTH 

PR-CM/KM  PR-CM  MB  PR-CM/KM  MB  PR-CM 

1 . 387E+00  2.845E+00  1.406E+03  2.354E+03 
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SPRINC-S 


(WX/W3-W1V)  UISN3Q  3QIX0IQ  N08W3 


Figure  29.  Comparison  of  Carbon  Dioxide  Density  for  the  Two 
Model  Atmosphere. 
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Figure  30. 


Comparison  of  Water 
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ALTITUDE  (KM) 


be  seen  in  Figure  45.  In  these  cases,  KjC  used  two  procedures. 

In  the  first  case,  there  were  no  clouds  above  1.5  km  in  the 
calculation,  and  in  the  second  case,  KSC  used  cloud  cover  at  3 
km.  The  effect  this  had  on  predicted  thermal  energy  at  a  given 
location  can  be  seer  in  Figures  46  thru  49. 

3.4  KSC  TRANSMISSION  TECHNIQUES. 

The  TAXV  code  was  developed  at  KSC  to  calculate  thermal 
environments  for  specific  cases  of  interest  using  the  Monte  Carlo 
data  base.  The  code  reads  in  a  specific  burst  altitude  and 
receiver  altitudes  and  ranges  of  interest  along  with  the  albedo 
surface  specification  and  the  surface  visibility  and  water  vapor 
density.  The  code  then  computes  the  irradiance  at  each  of  the 
receiver  points  as  discussed  previously  and  finds  the  exposure  by 
integrating  over  time. 

The  transmission  factor  given  previously  in  the  expression 
for  the  irradiance  is  given  by  the  product  of  the  direct  trans¬ 
mission  and  a  scattering  buildup 

T  -  VBs  ■  exP(-R*[  r  e<  X.H,)  -  Te<  A,  Hr)|  /  |nB-Hr  )*tir*bs  (31) 

where  the  exponential  factor  represents  the  attenuation  due  to 
aerosol  and  molecular  scattering  and  ozone  absorption,  TJR  refers 
to  the  molecular  absorption  in  the  IR  portion  of  the  spectrum, 
and  Bg  represents  the  buildup  due  to  multiple  scattering  in  the 
atmosphere.  The  quantity  re  is  an  optical  depth  defined  by 

50 

re(A,Hs)  -  fa  e (  A  ,H)  dH  (32) 
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where  oe  is  the  extinction  cross  section  (km-*)  given  by 


0e  (  X  *  H) 


oJX.H)  +  a  [  x  ,  h)  + 


aoz(*'H) 


(33) 


All  quantities  are  dependent  upon  wavelength  and  altitude  with 


V 


O  and 

K 


0QZ  being  the  aerosol,  Rayleigh  and  ozone  cross 


sections  respectively.  Ozone  is  a  molecular  absorber  and  the 
details  of  the  cross  section  will  be  given  in  a  later  section. 
The  ozone  absorption  varies  exponentially  with  ozone  concentra¬ 
tion  so  for  convenience  it  is  combined  with  the  scattering  cross 
sections  in  defining  the  model  atmosphere  parameters.  If  Hg  = 
Hr,  then  the  above  exponential  expression  is  replaced  by 


exp  (-R*  O  e  (X  ,hg) ) . 


(34) 


The  previous  two  sections  describe  how  the  direct  flux 
attenuation  is  computed  including  the  effects  of  the  scattering 
interactions  and  the  gaseous  absorption.  The  photons  scattered 
from  the  main  beam  can  undergo  multiple  scattering  interactions 
with  the  aerosol  and  molecular  constituents  of  the  atmosphere  and 
can  contribute  to  the  total  irradiance  at  the  receiver.  Various 
techniques  and  expressions  (Reference  11)  have  been  defined  in 
the  past  to  represent  this  scattered  buildup.  The  techniques 
based  upon  experimental  measurements  in  general  are  not  reliable 
when  extrapolated  beyond  the  range  of  the  original  experimental 
conditions. 


The  most  reliable  scattered  flux  predictions  are  obtained 
from  Monte  Carlo  Transport  calculations  (References  11  and  14) . 
Given  a  particular  meteorological  range  and  humidity  on  the 
ground  surface  then  the  scattered  flux  environment  can  be  deter¬ 
mined  as  a  function  of  receiver  range,  altitude  and  orientation 
for  a  particular  source  wavelength  and  altitude.  These  types  of 
runs  are  very  expensive  and  time  consuming  and  cannot  be  used 
directly  for  general  environment  predictions.  Instead,  the  data 
from  such  runs  must  be  processed  and  converted  to  a  form 
permitting  fast  accurate  environment  predictions. 

The  buildup  factor  approach  has  proven  to  be  the  best 
approach  for  this  type  of  prediction  and  has  been  used  in  many 
previous  transmission  predictive  techniques.  The  buildup  typi- 


R6 


cally  has  been  derived  from  elementary  transport  theory  or 
experimental  measurements  such  as  used  in  the  UK  approach.  The 
current  procedure  temporarily  uses  buildup  factors  taken  from  the 
RECIPE  transmission  routines  which  were  derived  to  match  the 
Monte  Carlo  data  bases  available.  The  buildup  factors  can  be 
defined  with  the  source  wavelength  and  altitude,  the  receiver 
range,  altitude  and  orientation  and  the  model  atmosphere  vari¬ 
ables  as  parameters. 

Previous  studies  (Reference  11)  have  shown  that  marked 
simplification  occurs  if  the  buildup  factors  are  expressed  as  a 
function  of  the  scattering  optical  depth  (SOD)  from  the  source  to 
the  receiver.  If  H  +  H  , 

S  i 

SOD  »  (0  A  (x  /Hs)  +  aR  (  X  ,Hs))  *  R,  (35) 

and,  if  H  »  H 
s  r 

SOD  *  R*  l(  rA(  X  ,Hs)-  rA(  x  ,Hr)  +  TR(X,Hr)-  H  X  ,Hr)||  Hs  -  Hr> 

where  all  variables  have  been  previously  defined.  In  earlier 
work  buildup  factors  as  a  function  of  SOD  defined  above  were 
compared  with  the  source  altitude,  receiver  altitude,  wavelength, 
and  various  meteorological  variables  parameters.  It  was  found 
that  the  wavelength  dependence  is  the  primary  parameter  of 
importance  in  representing  the  transport  results.  The  source  and 
receiver  altitudes  were  relatively  unimportant. 

The  wavelength  dependence  in  the  current  model  is  shown  in 
Figure  31  for  wavelengths  from  .32  to  3.8.  There  is  seen  to  be  a 
smooth  transition  in  the  buildup  to  smaller  values  as  the  wave¬ 
length  is  increased.  This  transition  follows  that  noted  in  the 
Monte  Carlo  data.  This  type  of  variation  is  expected  considering 
the  wavelength  dependence  of  the  scattering  functions  shown  in 
Table  12. 

The  scattered  flux  does  not  show  this  smooth  variation  with  wave¬ 
length  in  the  IR  region  because  of  the  extreme  wavelength 


(36) 
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Figure  31.  Comparison  of  Scattering  Build-up  Factors 
for  Various  Wavelengths. 
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dependence  of  the  gaseous  absorption.  Since  the  buildup  factor 
comes  from  the  ratio  of  the  total  to  direct  flux,  it  displays  the 
characteristics  of  the  scattering  and  not  the  absorption  inter¬ 
actions. 

The  curves  in  Figure  31  were  obtained  from  the  expressions 
used  in  the  updated  TAXV.  They  were  taken  from  the  RECIPE  code 
which  defines  the  buildup  factor  as  follows: 


Hr  <49  km 


BX  (  T  g)  -  1.3  *  exp 

(A  f  g)/z 

(37) 

where 

*S  “  *A  +  TR 

(38) 

2  ■  1  +  .3  exp  (-2 

rs) 

(39) 

and 

A  *  fl  exp  (-  y  V) 

(40) 

where  V  is  the  visual  range  (meteorological  range) .  The  variable 
is  a  function  of  the  field  of  view  (FOV)  given  by 

U  »  2.96  exp  (-2  FOV)  +  .055  exp  (-.05  FOV)  (41) 

and  a  is  a  function  of  wavelength  (m)  given  by 

a  -  . 38 X  ”'145  +  .025)  ”2,08.  (42) 

Hr  >  49  km 

BX  (  r  s)  -  B  exp  (A  T  s) 
where  A-  (.4  +  .1  +  2.29  v"*65  _1) 

-■■■■■  -  *  'i 

\(  1 . 0 8  +  .00127  V*,A) 

and  B  -  1  +  8  (.15/  X  )1,86exp  (.00926V)  (45) 


(43) 

(44) 
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It  is  assumed  that  the  orientation  dependence  of  the 
scattered  flux  is  the  same  as  the  direct  flux.  This  is  reason¬ 
able  because  the  scattering  is  strongly  forward  peaked  resulting 
in  a  narrow  angular  distribution  about  the  source  to  receiver 
direction. 

3.5  KSC  PROBABILISTIC  TRANSMISSION  METHOD. 

A  method  has  been  developed  to  compute  probabilistic  trans¬ 
mission  factors  and  exposures  using  the  KSC  thermal  atmospheric 
transmission  and  source  output  models  described  in  the  preceding 
sections.  This  methodology  is  detailed  in  Figure  32.  As 
discussed  previously,  the  first  step  in  the  procedure  is  a  care¬ 
ful  analysis  of  the  weather  at  the  region  of  interest.  The  data 
base  containing  the  build  up  factors  and  the  model  atmospheres 
stored  in  TAXV  must  be  compatible  with  the  altitude  dependent 
atmospheric  parameters  for  the  region  of  interest.  The  current 
models  as  discussed  previously  are  applicable  to  mid-latitude 
temporate  regions  with  a  fall-winter  or  spring-summer  seasonal 
choice.  Other  regions  such  as  the  tropic,  mideast  desert,  or 
arctic  regions  will  require  additional  modelling  of  atmospheres 
for  these  regions.  Necessary  information  is  available  in  the 
literature  and  in  fact  considerable  experience  for  the  arctic  and 
tropic  regions  has  been  obtained  in  previous  studies.  When  the 
data  base  is  expressed  as  buildup  factors  as  a  function  of 
scattering  optical  depth  relatively  small  variations  will  be 
noted  for  the  various  atmosphere  models.  Storing  the  data  as 
scattered  flux  as  a  function  of  the  geometrical  variables  as  is 
done  in  QRT,  results  in  a  data  bane  with  greater  variation 
between  model  types  and  with  no  well  defined  technique  for 
interpolating  results  between  models. 

The  weather  analysis  required  in  the  KSC  approach  involves 
defining  a  set  of  cases  involving  relevant  atmospheric  parameters 
each  of  which  represents  an  equal  occurrence  probability. 
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Typically  the  parameters  of  most  importance  are  the  visibility 
(or  meteorological  range)  and  the  albedo  surface  characteristics 
including  ground  elevation  and  albedo  and  cloud  ceiling  altitude 
and  coverage.  For  some  conditions  it  may  be  necessary  to  include 
additional  parameters  such  as  the  water  vapor  density.  These 
parameters  can  be  described  in  any  desired  detail  as  monthly 
seasonal  or  annual  averages. 

The  approach  used  in  the  UK  analysis  was  to  consider  cloud 
ceiling  altitude  and  visibility  as  parameters.  The  details  of 
the  UK  analysis  were  discussed  in  Section  2.  Basically,  the 
probability  of  occurrence  of  visibility  was  given  for  three 
different  ceiling  altitudes.  For  the  demonstration  calculations 
and  comparisons,  20  different  sets  of  weather  conditions  were 
defined  each  of  which  represents  an  occurrence  probability  of  5%. 
These  20  different  cases  are  shown  in  Table  16. 

A  serious  problem  with  the  UK  analysis  was  their  equating  a 
high  altitude  ceiling  with  no  cloud  cover.  These  imply 
completely  different  transmission  probabilities  as  will  be  shown 
in  the  following  discussion  and  led  to  the  definition  of  the  two 
alternate  sets  in  Table  16. 

There  is  a  certain  range  of  visibilities  associated  with 
'iach  condition.  The  maximum  visibility  for  each  condition  is 
used  to  identify  the  case.  Thus  this  gives  the  upper  limit  to 
the  transmission  associated  with  the  range  of  visibilities  for 
that  case.  This  is  the  proper  approach  for  USANCA  applications 
when  sure  safe  type  environments  are  being  developed.  The  i  and 
j  indices  are  merged  into  the  1  index. 

The  next  step  as  indicated  in  Figure  32  is  to  compute  the 
transmission  as  a  function  of  range  for  each  of  the  20  sets  of 
conditions.  This  is  straightforward  and  involves  running  the  TAXV 
code  for  each  case  and  storing  the  transmission  factor  as  a 
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WEATHER  ANALYSIS 

DEFINE  PROBABILITY  VS.  VISIBILITY 
CLOUD  LAYERS  DEPEND  ON  REGION 
MONTHLY,  SEASONAL,  ANNUAL 
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TRANSMISSION  CURVES  DIRECTLY  REPRESENT  OCCURRENCE  PROBABILITY 


Figure  32.  KSC  Probabilistic  Transmission  Method. 
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function  of  the  range.  One  then  has  a  matrix  of  transmissions  as 
shown  in  Table  17  for  the  20  different  cases  (1  index)  and  ranges 
(R^) .  This  table  represents  the  case  where  clouds  were  assumed 
to  exist  at  3km  for  group  III  conditions.  As  expected  the 
transmission  in  general  increases  as  the  visibility  increases 
within  each  albedo  group  at  a  particular  range.  A  priori  one 
cannot  predict  the  relative  transmissions  for  different  groups 
because  the  enhancement  factor  from  the  albedo  surfaces  depends 
upon  both  cloud,  source,  and  receiver  altitudes  and  the  range  in 
a  complicated  fashion.  In  general,  as  indicated  previously  for 
higher  ceilings,  the  maximum  enhancement  occurs  at  larger  ranges. 

These  results  are  plotted  in  Figure  33  where  the  symbols  A 
through  J  represent  cases  1  through  20,  respectively,  and  in 
Figure  34  where  the  data  points  are  connected  by  straight  lines. 
Recall  that  each  case  represents  a  51  occurrence  probability. 

The  next  step  as  noted  in  Figure  32  is  to  reorder  the  1 
index  independently  at  each  R^  in  order  of  increasing  trans¬ 
mission.  This  results  in  Figure  35  where  now  the  symbols  are 
associated  with  the  reordered  1  index  at  each  R^  and  represent 
increasing  (or  equal)  transmissions.  Since  each  case  represents 
an  occurrence  probability  of  5%,  the  reordered  1  index  at  each  R^ 
and  each  symbol  also  represents  an  integrated  probability  from 
.05  to  1  in  steps  of  .05.  The  data  are  shown  as  continuous 
curves  in  Figure  36  where  finally  the  cumulative  occurrence 
probability  transmission  factors  are  displayed.  As  indicated 
previously,  the  set  of  conditions  corresponded  to  the  choice  of  a 
cloud  ceiling  at  3km  altitude,  choosing  no  cloud  at  that  altitude 
(an  equivalent  condition  in  the  UK  weather  analysis)  results  in  a 
much  lower  transmission  as  will  be  noted  in  the  comparisons  given 
in  the  next  section. 
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Table  18.  Transmission  vs.  Range  and  Case*  Cookes  Weather  Data 
KSC  Transmission,  Clouds  at  3.0KM  (Group  III). 
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Figure  33.  Sample  Case.  Data  Identified  by  Weather  Conditions, 
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Figure  34.  Sample  Case.  Transmission  Curves  for  Weather 
Conditions. 
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Figure  35.  Sample  Case.  Data  Identified  by  Cumulative 
Occurrence  Probability. 
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Figure  36.  Sample  Case.  Transmission  Curves  for 
Occurrence  Probability. 
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SECTION  4 


COMPARISONS  OF  KSC  AND  UK  RESULTS 

In  this  section  comparisons  of  KSC  and  UK  results  are 
presented.  Early  in  the  study  it  was  recognized  that  the  UK 
method  disagreed  with  KSC  transmission  results.  The  UK  method 
could  not  be  recommended  for  general  USANCA  purposes  because 
crucial  parameters  of  importance  were  not  included.  However ,  the 
bulk  of  the  KSC  data  base  was  developed  prior  to  a  redefinition 
of  atmospheric  models  as  discussed  in  subsection  3.2.  A  more 
recent  data  base  developed  for  the  QRT  code  was  based  on  later 
model  atmospheres  but  needs  considerable  extension  and  manipula¬ 
tion  before  it  can  be  used  for  general  thermal  predictions  as 
will  be  discussed  in  Section  5. 

The  approach  used  in  the  current  study  was  to  generate 
transmission  factors  relevant  to  the  UK  weather  analysis  and 
pinpoint  the  important  parameters.  As  noted  in  the  last  section 
an  example  calculation  was  completed  and  probabilistic  trans¬ 
missions  were  calculated. 

Early  in  this  study  a  series  of  runs  were  made  with  TAXV 
using  parameters  obtained  from  the  UK  analysis.  The  code  at  that 
time  was  structured  with  a  table  of  buildup  and  enhancement 
factors.  The  various  parameters  varied  in  the  runs  are  noted  in 
Table  4.1  and  were  chosen  to  indicate  the  relative  importance  of 
the  parameters  for  USANCA  thermal  predictions. 

Two  different  yields  and  source  spectra  were  chosen  to  show 
the  effect  of  spectra  on  the  transmission  calculations.  Figures 
37  and  38  show  the  effect,  respectively,  of  yield  and  surface 
interaction  spectral  changes  on  the  trai  smission  for  a  surface 
burst  with  a  ground  albedo  with  nominal  dirt  values,  a  cloud 
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Table  19.  Transmission  Run  Parameters 


YIELDS:  10KT,  100KT 

GROUND  ALBEDO:  ZERO,  DIRT,  SNOW 

CLOUD  BASE  ALTITUDES:  NONE,  .3,  1.5,  3.0  KM 

SPECTRA:  SURFACE,  FREE  AIR 

SOURCE  ALTITUDE:  ZERO,  1KM 

VISIBILITY:  25KM,  6.5KM 

RECEIVER:  FLAT  PLATE  ON  GROUND  FACING  SOURCE 

HUMIDITY:  1.5,  10G/M* 
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Figure  37.  Effect  of  Yield  on  Transmission  Predictions. 
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Figure  38.  Effect  of  Spectrum  Differences  Caused  by 
Surface  Interactions  on  Transmission. 
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Range  (kia) 


ceiling  at  .3  km  and  a  visibility  of  25  km.  The  differences  are 
seen  to  be  negligible  in  both  cases.  Figure  39  shows  the  effects 
of  source  altitude  on  the  transmission.  The  transmission  for  the 
1  km  source  altitude  is  much  higher  than  for  surface  bursts. 

This  effect  is  not  allowed  for  in  the  UK  method  and  their 
transmission  corresponds  in  general  to  surface  bursts  and 
receiver  altitudes.  Neglect  of  this  dependence  can  seriously 
underpredict  the  thermal  exposure. 

Figures  40  through  43  show  the  effects  of  the  albedo 
surfaces  on  the  transmission  for  various  conditions.  In 
Figure  40,  the  transmission  factors  for  ceiling  altitudes  of  .3, 
1.5  and  3  km  are  compared  to  the  transmission  with  no  cloud  and 
with  zero  ground  albedo  assumed.  A  buildup  is  noted  in  each  case 
over  the  no  cloud  case.  The  range  at  which  the  maximum  buildup 
is  noted  increases  as  the  cloud  ceiling  is  raised.  Note  that  at 
large  ranges  the  transmission  with  a  cloud  ceiling  becomes 
smaller  than  for  the  no  cloud  case. 

Comparable  results  are  shown  in  Figure  41  for  a  dirt  ground 
surface  albedo.  The  differences  are  much  larger  because  of  the 
"channeling"  of  thermal  energy  to  much  longer  ranges  because  of 
the  reflective  surfaces.  In  Figure  42  the  effect  of  the  ground 
albedo  is  shown  for  the  .3  km  cloud  ceiling.  The  albedo  for  the 
dirt  surface  corresponds  to  sandy  soils. 

In  Figure  43,  the  combined  effects  of  the  albedo  surfaces 
and  burst  altitudes  are  shown  for  dirt  ground  albedo  and  25  km 
visibility.  Again  raising  the  burst  altitude  to  1  km  results  in 
a  dramatic  increase  in  the  transmission  factor. 

The  above  examples  all  used  a  25  km  visibility.  The  effect 
of  surface  albedo  for  a  6.5  km  visibility  is  noted  in  Figure  44 
and  should  be  compared  with  Figures  37-52  where  comparable 
results  are  given  for  a  25  km  visibility.  The  same  general 
trends  are  noted  but  as  expected  the  enhancement  is  less  for  the 
smaller  visibility  and  the  transmission  is  less  at  the  longer 
ranges. 


104 


Figure  39.  Effect  of  Source  Altitude  on  Transmission  Predictions. 
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Figure  40.  Effects  of  Cloud  Ceiling  Altitudes  on  the  Transmission 

for  Zero  Ground  Albedo. 
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Range  (km) 

Figure  41.  Effect  of  Cloud  Ceiling  Altitude  on  the  Transmission  for  Dirt 
Ground  Surface. 
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Figure  42.  Effect  of  Ground  Albedo  on  the  Transmission. 
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Figure  43.  Effect  of  Burst  Altitude  on  Transmission  with  Albedo 
Surfaces. 
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Figure  44.  Effect  of  Surface  Albedo  for  6.5  km  Visibility. 
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Range  (km) 


As  discussed  in  the  previous  section  a  set  of  probabilistic 
transmission  factors  was  developed  using  the  UK  weather  analysis. 
The  95%  and  the  50%  curves  are  compared  in  Figure  45.  The  solid 
lines  represent  the  KSC  analysis  with  no  clouds  above  1.5  km. 

The  dashed  lines  refer  to  the  KSC  analysis  assuming  clouds  at 
3  km  for  the  group  3  cloud  ceiling  condition  discussed  pre¬ 
viously.  As  expected  the  cases  with  clouds  have  a  much  higher 
transmission  than  the  cases  without  clouds.  Shown  also  are  the 
UK  curves  given  by  Cooke.  The  95%  curves  are  considerably 
different,  whereas,  the  Cooke  50%  curve  is  fairly  close  to  the 
KSC  curve  without  clouds  out  to  ranges  of  about  8  km.  The  Cocke 
95%  curve  does  not  appear  to  show  the  influence  of  mere  than  one 
cloud  layer  because  of  the  narrow  range  of  ground  ranges  where 
albedo  surface  enhancement  occurs.  The  Cooke  95%  curve  is  seen 
to  be  larger  than  the  KSC  curve  only  for  ranges  of  about  2  km  and 
for  ranges  of  10  km  and  greater.  The  flattening  out  of  the  Cooke 
95%  and  50%  curves  beyond  5  km  is  not  expected.  The  KSC  curves 
with  a  cloud  present  a  small  percentage  of  the  time  show  how 
crucial  the  difference  between  no  cloud  ceiling  and  a  cloud 
ceiling  at  3  km  is  even  though  both  cases  were  treated  as  equiva¬ 
lent  in  the  UK  analysis. 

Of  ultimate  concern  is  the  exposure  prediction  one  obtains 
for  particular  cases  of  interest.  In  Figure  46  the  thermal 
exposure  as  a  function  of  ground  range  is  shown  for  a  100  kt 
surface  burst  for  the  95%  case  with  the  curve  identification 
being  the  same  as  in  the  previous  figure.  All  curves  use  the 
same  source  output  as  described  in  Section  3.1  and  the  only 
difference  is  the  transmission  factors.  The  Cooke  transmission 
factor  leads  to  an  underprediction  of  the  exposure  for  exposure 
levels  less  than  about  30  cal/cma.  The  presence  of  clouds  at 
3  km  results  in  an  increase  in  the  exposure  for  levels  less  than 
about  9  cal/cma. 
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Figure  46.  Comparison  of  KSC  and  Cooke  Results  100  KT  Surface 
Burst:  95%  Case. 
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In  Figure  47  the  exposures  for  the  50%  case  are  shown.  The 
Cooke  curve  under  predicts  the  exposure  for  levels  less  than  8 
cal/cm2.  The  3  km  cloud  ceiling  predicts  higher  exposures  for 
all  levels  less  than  60  cal/cm2. 

The  same  exposure  predictions  are  shown  for  a  1  Mt  surface 
burst  in  Figures  48  and  49.  The  differences  are  emphasized 
relative  to  the  100  kt  case  and  begin  to  show  up  at  higher  expo¬ 
sures.  These  comparisons  indicate  that  the  Cooke  curves  should 
not  be  used  for  the  USANCA  applications.  The  KSC  curves  are 
uncertain  because  of  data  base  uncertainties.  As  indicated 
previously,  the  major  uncertainties  occur  at  small  visibilities. 
The  results  at  higher  visibilities  corresponding  to  larger 
probabilities  will  not  have  large  uncertainties.  In  the  next 
section,  the  QRT  data  base  is  described,  required  extensions  are 
discussed  and  sample  results  are  shown. 
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Figure  47.  Comparison  of  KSC  and  Cooke  Results  100  KT 
Surface  Burst:  50%  Case. 
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Figure  48. 


Comparison  of  KSC  and  COOKE  Results  1  MT  Surface 
Burst:  95%  Case. 
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Figure  49.  Comparison  of  KSC  and  COOKE  Results  1  MT  Surface 
Burst:  50%  Case. 
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SECTION  5 


MONTE  CARLO  DATA  BASE  DEVELOPMENT 


As  discussed  the  QRT  data  base  was  generated  using  recent 
model  atmospheres  although  only  one  aerosol  distribution  (rural) 
was  « ctually  used  in  defining  aerosal  scattering  phase  functions. 
The  data  base  was  generated  with  a  fairly  limited  set  of  meteoro¬ 
logical  parameters  and  other  limitations  which  restrict  the 
usefulness  without  considerable  modification  and  extension.  In 
this  section,  the  data  base  is  described,  and  the  limitations  are 
identified. 

5.1  QRT  DATA  BASE  AND  PREDICTIVE  APPROACH. 

In  Table  20,  the  major  parameters  of  the  QRT  data  base  are 
identified.  The  data  base  basically  contains  3  groups  of  data. 
Data  sets  1,  2,  and  3  contain  data  for  the  3  indicated  meteorolo¬ 
gical  ranges  for  two  different  ground  albedos  and  no  clouds. 

Data  sets  4  and  5  contain  data  for  the  same  geometry  except  the 
ground  altitude  is  at  1.524  km  (5  kft) .  Data  sets  6  and  7 
contain  data  for  two  different  cloud  ceiling  altitudes  as  a 
function  of  a  ground  albedo.  All  the  data  are  raw  unsmoothed 
Monte  Carlo  data  representing  the  scattered  flux  received  on  a 
flat  plate  detector  oriented  at  119  ((15-2)*  9+2)  different 
possible  directions.  Only  six  different  wavelengths  were  used 
which,  as  noted  previously  (and  detailed  in  the  following  sub¬ 
section)  ,  can  not  adequately  represent  the  nuclear  thermal 
spectrum. 

The  QRT  approach  consists  of  calculating  the  direct  flux  at 
six  wavelengths  including  scattering  and  narrow  bank  IR  absorp¬ 
tion  effects.  The  scattered  flux  is  determined  by  interpolating 
in  the  massive  Monte  Carlo  data  base  which  stores  the  scattered 
flux  per  unit  source  versus  source  altitude,  receiver  altitude 
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Table  20.  QRT  Data  Base 
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and  range,  azimuth  and  polar  angles,  ground  albedo,  and  ground 
level  meteorological  range.  The  fluxes  are  curve  fit  at  up  to  20 
wavelength  intervals.  The  receiver  exposure  is  then  calculated 
by  integrating  the  direct  and  scattered  contributions  over  the 
fireball  spectral  output. 

The  direct  flux  from  a  point  source  is  represented  by 


(X,) 


008  -r  / cosff  * 
TZtr-  e  °  fa 


(X,  X) 


where : 


(46) 


x 

*z 
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Fa  (XfX) 


source  to  receiver  distance 

angle  between  the  source  to  receiver  vector  and 
the  receiver  surface  normal 

angle  between  the  source  to  receiver  vector  and 
the  ground  surface  normal 

atmospheric  scattering  optical  thickness  between 
the  source  and  receiver  planes 
gaseous  attenuation  factor  for  wavelength  and 
source-receiver  distance  x. 


If  the  source  and  receiver  are  at  the  same  altitude  then  the 
exponential  argument  is  replaced  by  -sx  where  s  is  the  atmo¬ 
spheric  scattering  extinction  coefficient  at  the  source  altitude. 

As  described  in  the  QRT  manual,  the  gaseous  attenuation 
factor  accounts  for  the  attenuation  by  atmospheric  gases  in  the 
infrared  wavelength  region.  Gases  included  in  the  modeling 
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include  H20,  CO^,  CH4 ,  and  0 2  for  the  wavelengths  chosen  in  the 
data  base.  The  transmission  along  a  path  is  given  by 


fa  -  exP 


(AH  (j,  hx)  -  AH(j,  h2))/|cos0|  (47) 

j  t/f"+  ^AH(j,  h1}  -  AH(j,h2)  )  */  icosgl 

><y  (BH ( j ,  hx)  “  BH  ( j  ,  h2)) 


2  0  ( 


where  h^  and  h2  are  the  source  and  receiver  altitudes  respec¬ 
tively,  is  defined  above,  and  the  summation  is  over  the 
contributing  absorbing  species  for  the  various  wavelengths  (H20 
for  .94  pm,  H20,  C02,  and  02  for  1.25  pm  and  H20,  C02,  and  CH4 
for  1.7  pm).  FA  =  1  for  wavelengths  of  .337,  .55  and  1.06  pm. 
The  AH  and  BH  arrays  are  defined  by 

CD 

AH ( j ,  h)  =  J  A ( j ,  h)  *  m( j ,  h) dh  (48) 

•h 


and 


BH ( j ,  h) 


_ 4 _ 

1013  * 


00 

J  BMj,  h)  *  m(h,  h)  *  p (h) dh,  (49) 

h 


where  m  is  the  absorber  density  and  p  is  the  atmospheric  pres¬ 
sure.  A  and  B  are  parameters  defined  in  the  Goody  model  and  are 
given  by 


A 


and 


(50) 


where  S(i)  is  the  line  strength  and  a(i)  is  the  line  width  and 
the  summation  is  over  the  lines  in  the  relevant  wavelength 
interval. 


This  approach  is  based  upon  a  wavelength  interval  of  200 
wavenumbers  which  is  the  basis  interval  used  in  the  AFGL 
absorption  data  base.  The  actual  width  of  the  IR  intervals  is 
considerably  larger  than  200  (as  will  be  discussed  later)  which 
can  lead  to  inaccuracies  in  the  transmission  predictions  in  the 
IR  region. 
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Note  that  0^  was  not  included  as  an  absorbing  species  in  the 
IR  and  apparently  was  not  included  in  the  UV  and  visible  portion 
of  the  spectrum.  0^  can  be  important  in  that  portion  of  the 
spectrum  and  can  be  easily  included  in  the  a  and  T  variables 
in  the  exponential  argument  in  the  direct  flux  expression. 

The  scattered  flux  contributions  are  calculated  by  inter¬ 
polating  in  the  massive  data  base.  The  data  represent  the  energy 
incident  upon  a  small  plane  receiver  which  is  oriented  in  parti¬ 
cular  directions  at  a  particular  location  in  space.  Data  are 
included  for  specified  source  altitudes  and  for  certain 
meteorological  ranges  and  albedo  surface  conditions  for  the  six 
defined  wavelengths.  The  data  are  the  results  of  a  large  number 
of  Monte  Carlo  runs  for  the  particular  conditions  of  interest  in 
which  angular  sampling  was  performed  to  define  the  scattered 
intensity  as  a  function  of  angle.  The  scattered  intensity  was 
then  converted  to  the  more  useful  results  stored  in  the  data 
base. 


The  data  base  parameters  were  summarized  in  Table  20.  Seven 
different  data  sets  are  available  for  the  particular  albedo 
surface  combinations  and  meteorological  ranges  shown.  For  each 
data  set  (1  to  5)  two  basic  scattered  flux  arrays  corresponding 
to  zero  ground  albedo,  and  the  additional  scattered  flux  for  an 
albedo  of  .8  are  stored  for  all  possible  combinations  of  the 
other  parameters  listed.  Each  data  set  contains  30  different 
data  groups  containing  the  results  of  each  of  the  six  wavelengths 
run  at  the  five  different  source  altitudes.  For  each  group 
scattered  data  are  stored  for  70  different  receiver  locations 
corresponding  to  the  seven  receiver  altitudes  and  the  ten 
receiver  ranges.  At  each  location  the  scattered  flux  is  stored 
as  a  function  of  thirteen  polar  and  nine  azimuthal  directions 
with  the  nadir  and  zenith  polar  directions  being  stored 
separately  since  the  azimuthal  directions  are  degenerate  for 
those  directions.  As  mentioned  arrays  are  defined  for  zero  and 
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non-zero  albedo  surfaces.  The  scattered  flux  for  a  unique 
combination  of  these  different  parameters  is  found  by  a  compli¬ 
cated  search  and  interpolation  scheme  for  each  of  the  wave¬ 
lengths.  If  the  particular  receiver  altitude  and  range  do  not 
correspond  to  the  values  in  the  table,  values  for  the  exposure 
are  computed  for  each  of  the  four  locations  braketing  the  loca¬ 
tion  of  interest  and  interpolation  is  then  used  to  define  the 
exposure  for  the  location  of  interest.  Interpolation  using  the 
cosines  of  the  polar  and  azimuthal  angles  is  performed  to  deter¬ 
mine  the  scattered  flux  valuec. 

As  discussed  previously  arrays  are  stored  representing  the 
value  of  the  scattered  flux  for  an  albedo  value  of  zero  and  the 
difference  resulting  from  an  albedo  of  0.8.  An  approximate 
procedure  is  then  given  to  compute  the  scattered  flux  for  any 
arbitrary  albedo. 

Using  the  methods  described  above,  values  of  the  direct  and 
scattered  flux  are  computed  for  each  of  the  six  wavelengths.  The 
albedo  correction  is  done  separately  for  each  wavelength.  The 
direct  flux  and  the  scattered  flux  are  then  fit  separately  as  a 
function  of  wavelength.  As  indicated  in  Figure  50,  the 
0.337  /im,  0.55  urn  and  the  1.06  #*m  results  are  fit  to  evaluate 

fluxes  for  wavelengths  less  than  1.06  with  the  exception  of  the 
region  from  .926  to  .978  where  the  .94  results  are  treated 
separately.  The  narrow  bands  between  1.1  and  1.15  |im  and  1.35 
and  1.5  #im  represent  strong  absorption  bands  and  the  energy  in 
these  regions  are  set  to  zero.  The  1.25  and  1.7  fi  m  results  are 
used  throughout  their  respective  bands  as  shown.  This  curve 
illustrates  one  of  the  penalties  of  using  the  scattered  flux 
approach.  *  In  the  IR  region  the  scattered  flux  as  well  as  the 
direct  flux  shows  a  wide  and  non-predictable  dependence  upon 
wavelength  because  of  the  complicated  IR  absorption  characteris¬ 
tics.  The  curve  fitting  procedure  cannot  be  used  consistently 
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Wavelength  (  m) 


across  the  entire  wavelength  region  and  the  arbitrary  approach 
described  above  must  be  used.  The  buildup  approach  described 
previously  which  has  been  used  in  all  of  the  KSC  thermal  trans¬ 
port  calculations  is  a  much  more  reasonable  and  consistent 
approach  and  can  be  used  with  confidence  over  the  entire  wave¬ 
length  region. 

Given  the  curve  fits  as  described  above,  the  fluxes  are 
evaluated  at  up  to  20  wavelengths  and  folded  in  with  a  time 
dependent  spectral  source  output.  Integration  over  wavelength 
then  gives  the  irradiance  at  the  receiver  and  integration  over 
time  gives  the  exposure. 

A  very  large  data  base  is  thus  available  in  the  QRT  code  but 
it  has  been  found  that  significant  limitations  in  the  predictive 
approach  and  data  base  exist  and  that  significant  modifications 
are  required. 

5.2  LIMITATIONS  IN  THE  QRT  CALCULAT I ONAL  APPROACH  AND  DATA  BASE. 

5.2.1  Inadequate  Wavelength  Base. 

As  discussed  in  Section  5.1,  there  should  be  additional 
wavelengths  in  the  IR  region  to  represent  properly  the  nuclear 
weapon  thermal  output  spectrum  adequately  in  this  region.  Adding 
these  wavelengths  and  converting  to  the  buildup  factor  approach 
to  calculate  the  scattered  flux  would  provide  the  best  approach 
to  solve  this  problem. 

5.2.2  Utilization  of  Scattered  Flux  Data. 

As  discussed  previously,  the  scattered  flux  magnitude  varies 
very  rapidly  with  wavelength  in  the  IR  region  because  of  the  very 
complicated  molecular  absorption  band  structure  of  the  absorbing 
gases.  In  absorbing  regions,  the  scattered  flux  will  be  very 
small  relative  to  that  in  the  windows  adjacent  to  the  absorbing 


125 


bands.  Trying  to  fit  the  scattered  flux  as  a  function  of 
wavelength  is  possible  only  in  the  UV  and  visible  portion  of  the 
wavelength  region  as  was  shown  in  the  previous  section.  In  the 
IR  region,  it  was  found  necessary  to  use  band  averaged  values. 

To  represent  accurately  the  wavelength  dependence  of  the 
scattered  flux  would  require  many  more  source  wavelengths  than 
currently  used.  This  problem  is  solved  by  using  buildup  factors 
as  described  in  subsection  3.2. 

5.2.3  Utilization  of  Unsmoothed  Monte  Carlo  Data. 

The  scattered  flux  data  stored  on  the  tapes  are  basically 
unsmoothed  Monte  Carlo  data  that  have  been  integrated  over 
sampling  angles  to  represent  the  flux  on  flat  plate  receivers 
with  a  particular  orientation  at  a  particular  location  in  space. 
Monte  Carlo  data  always  have  statistical  uncertainties  which 
makes  their  usage  suspect  without  using  smoothing  or  curve 
fitting  techniques.  Large  variations  can  occur  in  the  scattered 
data  as  a  function  of  receiver  range  and  altitude  and  orientation 
angle.  These  variations  can  lead  to  inconsistent  trends  of  the 
exposure  versus  range.  Interpolation  in  the  very  large  data  base 
is  extremely  time  consuming.  Curve  fitting  and  smoothing  the 
buildup  factors  is  a  much  more  efficient  and  safe  approach. 

5.2.4  Narrow  Band  IR  Absorption  Parameters. 

IR  absorption  parameters  averaged  over  a  band  width  of 
200  wavenumbers  were  used  in  the  model.  In  Table  21,  the  width 
of  the  bands  as  discussed  in  QRT  supporting  documentation  is 
shown.  The  .337  and  the  .55  wavelengths  are  not  treated  as  bands 
in  the  modeling.  All  of  the  bands  are  seen  to  be  much  wider  than 
the  200  wavenumbers  used  to  define  the  basic  absorption  para¬ 
meters.  Significant  errors  exist  in  representing  the  wide  band 
absorptions  by  narrow  bands  at  the  band  center. 
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Table  21.  QRT  Band  Structure 


Band 

Band 

Band 

Wavelength 

Limits 

Width  ■ 

(  *-m) 

( 

(wavenumber) 

.337 

N/A 

N/A 

.55 

N/A 

N/A 

.94 

.926  -  .978 

574 

1.06 

.978  -  1.10 

1134 

1.25 

1.15  -  1.35 

1288 

1.70 

1.50  -  1.80 

1111 
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SECTION  6 


RECOMMENDED  APPROACH  FOR  QRT  DATA  BASE  REVISIONS 

In  this  section  the  recommended  approaches  for  using  the  QRT 
code  data  base  are  discussed.  The  KSC  wavelength  band  structure 
has  been  found  to  represent  the  output  spectrum  and  if  the 
buildup  factor  approach  is  used  will  give  a  very  good  description 
of  the  atmospheric  transport  properties.  This  basic  wavelength 
band  structure  will  be  used  with  revised  IR  absorption  techniques 
to  calculate  the  direct  flux  as  described  in  detail  in 
Section  3.4. 

The  irradiance  as  a  function  of  time  will  be  given  by  the 
expression 

H  =  g  Y  T  (X,  h  ,  h  ,  R)  *  B  (A,  h  ,  h  ,p  ,*)  (51) 

4 1  Rs  x  • 

*  B  (X,  h  ,  h  ,  R,  #)  *  P(  A  ,  r,  cos0) 

Q 

where  TD  =  the  direct  flux  transmission  discussed 

previously 

Bg  =  buildup  factor  with  zero  albedo 

Ba  =  buildup  enhancement  factor  due  to  albedo 
surfaces 

p  =  power  output  of  the  fireball 

A  =  wavelength 

h  ,  h  =  source  and  receiver  altitudes 

s '  r 

g  =  geometry  factor 

R  =  slant  range 

p  =  optical  depth  along  slant  path 

t  =  time  after  burst 

0  =  pol  ir  angle  to  slant  path 

0  =  orientation  of  sampling  surface 


128 


The  source  altitude  at  zero  time  is  the  height  of  burst  and 
thereafter  is  the  fireball  altitude  which  increases  as  a  function 
of  time.  The  receiver  altitude  and  the  slant  range  can  also  be  a 
function  of  time  for  a  moving  target. 

The  Bg  ( X  ,  hg,  hr>  R,  p  ,  0  )  represents  the  buildup  factor 

evaluated  for  all  of  the  listed  parameters.  These  will  be  stored 

as  curve  fits  in  as  many  parameters  as  possible  in  order  to 

reduce  the  running  time  and  interpolation  procedures  will  be  used 

in  the  other  parameters.  The  B  (  X  >  h  ,  h  ,  R,  0  )  represent 

h  s  r 

the  buildup  enhancement  factors  due  to  the  albedo  surfaces  and 
again  will  be  curve  fit  in  as  many  parameters  as  possible.  The 
P  (X  ,  t,  cos  0  )  represents  the  power  output  of  the  fireball  and 
will  be  obtained  from  the  revised  RECIPE  code  described  in 
Section  3.1.  The  exposure  is  then  found  by  integrating  the 
irradiance  over  time. 

6.1  EXTENSION  OF  THE  WAVELENGTH ’ DATA  BASE. 

As  discussed  in  Section  3.1  the  current  six  wavelengths  do 
not  extend  to  high  enough  wavelengths  to  represent  adequately  the 
nuclear  weapons  termal  output  spectrum.  In  order  to  calculate 
the  buildup  factors  for  the  entire  range  of  wavelengths  in  the 
KSC  wavelength  band  structure  needed  in  the  expression  given 
above  it  will  be  necessary  to  generate  scattered  flux  data  at 
longer  wavelengths  than  currently  are  in  the  QRT  data  base.  In 
order  to  make  maximum  usage  of  the  current  data  base  it  is 
suggested  that  two  additional  wavelengths  be  defined  at  2.7  and 
3.8  /zm.  Because  of  the  predictable  variation  of  the  buildup, 
factors  with  wavelengths  as  shown  in  the  previous  section,  it 
will  be  possible  to  represent  accurately  the  multiple  scattering 
buildup  with  the  resulting  eight  wavelengths.  Additional  Monte 
Carlo  runs  will  be  made  for  those  two  wavelengths  with  the  same 
general  geometrical  parameters  as  the  current  QRT  data  base  for 
the  .01,  5.0  and  80  kft  source  altitudes.  If  consistent  predict- 
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able  trends  are  noted,  it  may  not  be  necessary  to  make  runs  at 
all  of  the  possible  combinations. 


6.2  CONVERT  SCATTERED  FLUX  DATA  BASE  TO  BUILDUP  FACTORS. 


The  buildup  factors  will  be  generated  from  the  current  data 
base  by  dividing  the  scattered  flux  arrays  (SINT,  SIT,  and  SIB) 
at  each  of  the  receiver  locations  by  the  proper  element  of  the 
direct  flux  array  RES.  For  each  of  the  receiver  locations  the 
scattering  optical  depth  can  be  calculated  by  the  following 
equation 


SOD  ■ 

where  t  = 


and 


(52) 


the  integral  of  the  optical  depth  to  the 
altitude  of  interest 
source  and  receiver  altitudes 
the  slant  range. 


Example  plots  of  the  buildup  factor  for  a  wavelength  of  .55  fi m,  a 
source  altitude  of  .01  kft,  and  a  meteorological  range  of  23  km 
obtained  from  data  set  1  are  shown  in  Figure  51.  The  results  are 
shown  for  receiver  altitudes  of  .001,  1.0,  2.0,  and  5.0  km. 

These  factors  were  derived  for  a  receiver  orientation  facing  the 
source  by  interpolation  in  the  angular  information  in  the  data 
base.  Note  that  the  buildup  factors  for  the  various  altitudes  do 
not  differ  by  a  large  amount  when  plotted  as  a  function  of  the 
optical  depth.  The  points  plotted  were  calculated  from  the  ter. 
sampling  ranges  used  in  the  data  base.  For  the  lower  altitudes, 
a  larger  optical  depth  is  traversed  along  the  slant  path.  The 
maximum  optical  depth  plotted  is  seen  to  decrease  as  the  altitude 
is  increased.  If  the  horizontal  range  were  used  as  the 
independent  variable  the  curves  would  be  separate  and  distinct. 
Likewise  if  the  scattered  flux  were  treated  as  the  primary 
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Figure  51.  Plot  of  Buildup  Factor  with  Source  Altitude  .01  KFT. 
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dependent  variable  separate  curves  would  be  obtained  since  the 
flux  would  be  reduced  much  more  for  the  lower  altitudes  than  for 
the  higher  ones. 

Similar  results  are  plotted  in  Figure  52  for  the  2  kft 
source  altitude.  Again  note  that  the  curves  for  the  various 
sampling  altitudes  lie  close  to  the  same  line.  Comparing  the  two 
figures  one  notes  also  that  the  results  for  the  two  different 
source  altitudes  are  very  close  together  when  plotted  in  this 
manner.  This  indicates  that  representing  the  results  in  this 
manner  is  a  natural  approach  as  one  expects  from  elementary 
transport  theory. 

A  similar  procedure  would  be  used  in  deriving  the  enhance¬ 
ment  factors  from  the  albedo  surfaces  using  the  TINT,  RIT,  and 
RIB  arrays.  Previous  studies  at  KSC  using  the  TRAX  results  have 
indicated  that  in  this  case  the  more  natural  procedure  is  to 
present  the  results  as  the  buildup  enhancement  as  a  function  of 
the  distance  from  one  of  the  albedo  surfaces  and  the  horizontal 
range  with  the  unit  of  distance  being  the  distance  between  the 
surfaces.  This  indicates  that  the  effect  is  more  nearly  depen¬ 
dent  upon  the  geometry  of  the  problem  and  not  on  the  optical 
depths  involved. 

6.3  SMOOTHING  AND  CURVE  FITTING  BUILDUP  FACTORS. 

All  Monte  Carlo  results  are  statistically  uncertain  with  the 
variance  upon  the  penetration  depths  involved  in  the  problem,  tne 
efficiency  of  variance  reduction  techniques  used  in  the  code,  and 
the  number  of  histories  used  in  generating  the  results  which 
directly  impacts  the  running  times  of  the  prob’ems.  In  practice 
one  is  always  dealing  with  significant  statistical  uncertainties 
which  if  not  handled  carefully  can  lead  to  inconsistent  results. 
For  instance,  observe  the  pluses  plotted  in  Figure  52  repre¬ 
senting  a  receiver  altitude  of  5  km.  An  increase  in  the  buildup 
is  noted  as  the  SOD  is  increased  as  is  expected  from  transport 
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Figure  52.  Plot  of  Buildup  Factors  with  Source  Altitude  2KFT 


133 


SOD 


theory.  However,  at  about  .5  the  reverse  trend  is  noted  for 
three  points,  then  the  remaining  two  points  increase  with 
increase  in  SOD.  This  trend  is  not  real  but  merely  indicates  the 
uncertainty  in  the  results.  If  the  raw  Monte  Carlo  results  were 
used  with  interpolation  then  a  fallacious  dependence  of  the  flux 
versus  range  would  be  obtained.  Likewise  the  circles  show  a 
strange  unnatural  behavior  greater  than  1.5.  In  order  to  obtain 
consistent  results  a  cross  plotting  and  smoothing  procedure 
should  be  used. 

Even  if  the  data  are  smoothed  and  reduced  to  a  common  mesh 
table  lookup  and  interpolation  techniques  are  time  consuming.  A 
better  approach  is  to  curve  fit  the  data  in  as  many  parameters  as 
possible.  This  procedure  smoothed  results  and  is  much  faster 
than  the  other  approach.  As  discussed  in  Section  3.4,  curve  fits 
have  been  developed  which  are  capable  of  fitting  buildup  factors 
as  a  function  of  the  various  parameters  and  will  be  used  in  the 
future. 

6.4  BROAD  BAND  IR  ABSORPTION. 


Narrow  band  absorption  expressions  should  not  be  used  to 
represent  broad  band  transport  calculations  as  is  currently  done 
in  QRT  according  to  the  supporting  documentation  without  a 
thorough  investigation  of  the  influence  of  the  band  width  which 
has  not  been  done.  The  expressions  used  in  TRAX  were  derived 
from  experiments  in  the  laboratory  and  in  the  atmosphere.  The 
narrow  band  200nm  AFGL  expressions  discussed  earlier  have  been 
integrated  over  the  broad  bands  and  compared  with  the  current 
TRAX  expressions.  No  large  discrepancies  were  noted;  however,  in 
order  to  make  the  predictive  codes  consistent  with  AFGL  codes, 
such  as  low  TRAX,  the  integrated  AFGL  results  should  be  used. 
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